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FEDERAL  AVIATION  ADMINISTRATION 
SYSTEMS  RESEARCH  AND  DEVELOPMENT  SERVICE 
SPECTRUM  MANAGEMENT  STAFF 

Statement  of  Mission 

f 

The  mission  of  the  Spectrum  Management  Staff  is  to  assist  the  De- 
partment of  State,  Naitioriai  Telecommunications  and  Information 
Administration,  and  the  Federal  Communications  Commission  in  as- 
suring the  FAA's  and  the  nation's  aviation  interests  with  sufficient 
protected  electromagnetic  telecommunications  resources  throughout 
the  world  to  provide  for  the  safe  conduct  of  aeronautical  flight  by 
fostering  effective  and  efficient  use  of  a natural  resource- -the 
electromagnetic  radio  frequency  spectrum.  . . - - , 

This  object  is  achieved  through  the  following'servrcTesT 

. Planning  and  defending  the  acquisition  and  retention 
of  sufficient  radio  frequency  spectrum  to  support  the 
aeronautical. interests  of  the  nation,  at  home  and  a- 
brb'ad,  V and. spectrum  standardization  fop,  the;  . World's 
aviation  community i -•  * ’ i 

. Providing  research,  analysis,  engineering,  and  evalu- 
ation in  the  development  of  spectrum  related  policy, 
planning,  standards,  criteria,  measurement  equipment, 
and  measurement  techniques. 

. Conducting  electromagnetic  compatibility  analyses  to 
determine  intra/inter-system  viability  and  design 
parameters,  to  assure  certification  of  adequate  spec- 
trum to  support  system  operational  use  and  projected 
growth  patterns,  to  defend  aeronautical  services  spec- 
trum from  encroachment  by  others , and  to  provide  for 
the  efficient  use  of  the  aeronautical  spectnum. 


Developing  automated  frequency  selection  computer  , 
programs/routines  to  provide  frequency  planning,  fre- 
quency  assignment,  and  spectrum  analysis  capabilities 
in  the  spectrum  supporting  the  National  Airspace  Sys- 
tern.  Oj 

Providing  spectrum  management  consultation,  as,sis“. 
tance,  and  guidance  to  all  aviation  inferest,  users, 
and  providers  of  equipment  and  services,,  both  na- 
tional and  international. 
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PROPAGATION  MODEL  (0.1  to  20  GHz)  EXTENSIONS 
FOR  1977  COMPUTER  PROGRAMS 
G.  D.  GIERHART  and  M.  E.  JOHNSON ^ 

1.  INTRODUCTION 

Assignments  for  aeronautical  radio  in  the  radio  frequency 
spectrum  must  be  made  so  as  to  provide  reliable  services  for  an 
increasing  air  traffic  density  [19] Potential  interference  be- 
tween facilities  operating  on  the  same  or  on  adjacent  channels 
must  be  considered  in  expanding  present  services  to  meet  future 
demands.  Service  quality  depends  on  many  factors,  including  the 
desired-to-undesired  signal  ratio  at  the  receiver.  This  ratio 
varies  with  receiver  location  and  time  even  when  other  parameters, 
such  as  antenna  gain  and  radiated  powers,  are  fixed. 

In  1973,  an  air/ground  propagation  model  developed  at  the 
Department  of  Commerce  Boulder  Laboratories  (DOC-BL)  by  the  Insti- 
tute for  Telecommunication  Sciences  (ITS)  for  the  Federal  Aviation 
Administration  (FAA)  was  documented  in  detail.  This  IF- 73  (ITS- 
FAA-1973)  propagation  model  has  evolved  into  the  IF-77  model, 
which  is  applicable  to  air/air,  air/ground,  air/satellite,  ground/ 
ground,  and  ground/satellite  paths.  The  IF-77  has  been  incorpor- 
ated ir.to  10  computer  programs  that  are  useful  in  estimating  the 
service  coverage  of  radio  systems  operating  in  the  frequency  band 
from  0.1  to  20  GHz.  These  programs  may  be  used  to  obtain  a wide 
variety  of  computer-generated  microfilm  plots.  A plotting  capa- 
bility summary  is  provided  in  table  1,  and  program  input  para- 
meters are  summarized  in  tables  2 through  4.  These  tables  were 

^The  authors  are  with  the  Institute  for  Telecommunication  Sciences, 
National  Telecommunications  and  Information  Administration,  U.  S. 
Department  of  Commerce,  Boulder,  Colorado  d0303. 

^References  are  listed  alphabetically  by  author  at  the  end  of  the  ' 
report  so  that  reference  numbers  do  not  appear  sequentially  in 
the  text. 
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Table  1. 

Capability 


lolling** 

6 

UIBINS 

licflection  coefficient** 

7 

LOBING 

Path  length  difference** 

8 

LOBING 

Time  l.ig** 

9 

LOBING 

Lohtng  frcqiicncy-I)** 

10 

LOBING 

lol'ing  frcquency-H** 

11 

LOBING 

Reflrction  point** 

12 

LOBING 

Llcvation  tmglc** 

13 

LOBING 

Flcvation  angle  difference** 

14 

LOBING 

.Spectral  plot** 

IS 

LOBING 

I'over  available 

16 

ATCIA 

Power  dciuiity 

17-19 

ATOA 

Tran.'onission  loss 

20 

ATTIA 

Power  av.i liable  curves 

21 

ATLAS 

l^iwor  density  curves 

22 

ATLAS 

Iransmission  loss  curves 

23 

ATLAS 

Power  available  vol'jae 

24 

IIIKX) 

Power  density  volux- 

’S 

IIIKX) 

Iransiis.sion  loss  voluae 

26 

IIIPOI) 

1 IKP  coiitoiii's 

27-29 

AFOOS 

IViwer  available  contouirs 

30 

APOOS 

I'owvr  dciKity  contours 

31 

i\P0C6 

fransmission  loss  contours 

32 

•\P00S 

Signal  ratio-S 

3.T 

ATAai 

Plotting  Capability  Guide  [21,  table  1] 

Fl£ur«(a)*  Pnitf  . '*  Wwaato 

6 lOBtNS  TrwMtissian  loss  vorsus  path  distance. 

7 LOBINC  Effactiva  specular  reflection  coefficient  versus  path 

distance. 

8 LOBING  Difference  in  xaflaetad  «d  dixact  ray  lengths  versus 

path  distance. 

9 LOBING  Sane  as  above  with  path  length  difference  expressed  as 

tine  delay. 

10  LOBING  Nomelized  distance  lobing  frequency  versus  path  dis- 

tance. 

11  LOBING  Norsalized  height  lobing  frequency  versus  path  distance. 

12  LOBING  Distance  to  reflection  point  versus  path  distance. 

13  LOSING  Direct  ray  elevation  angle  versus  path  distance. 

14  LOBING  Angle  by  whieJ:  the  direct  ray  exceeds  the  reflected  ray 

versus  path  distance. 

15  LOBING  A^>litude  versus  frequency  response  curves  for  various 

path  distances. 

16  AIQA  Power  available  at  receiving  antenna  versus  path  dis- 

tance or  central  angle  for  tine  availirijilities  S,  SO, 
and  9S  percent. 

17-19  ATQA  Sinilar  to  above,  but  with  power  density  ordinate. 

20  ATOA  Similar  to  above,  but  with  transnission  loss  ordinate. 

21  ATLAS  Power  available  curves  versus  distance  are,  provided 

for  several  aircraft  altitudes  with  a selected  tine 
availability,  and  a fixed  lower  antenta  height. 

22  ATLAS  Sinilar  to  above,  but  with  power  density  as  ordinate. 

23  ATLAS  Sinilar  to  above,  but  with  transnission  loss  as  ordinate. 

24  IIIKX)  Fixed  power  available  contours  in  the  altitude  versus 

distance  plane  for  time  availabilities  of  5,  5b,  and 
95  percent. 

IIIPOI)  Sinilar  to  above,  but  with  fixed  power  density  contours. 

20  IIIPOI)  Similar  to  above,  but  with  fixed  transnission  loss 

contours. 

27' 29  APOOS  Contours  for  several  HIRP  levels  needed  to  ncet  a par- 

ticular power  density  requircannt  are  shown  in  the  ai- 
titufe  versus  distance  plane  for  a single  tine  availa- 
bility. 

3U  APOOS  Sinilar  to  above,  but  with  power  available  contours  for 

a single  EIRP. 

31  j\P006  Similar  to  above,  but  with  power  density  contours. 

32  iTPOns  Sinilar  to  above,  but  with  trannaission  loss  contours. 

■V  ATAHI  Uesiied-to-isidesired,  DA),  signal  ratio  vei^  station 

separation  for  a fixed  desired  facility-to-rcceiver 
dist.-uKCi  and  tiRK  availabilities  of  5,  50,  and  95 
percent. 


Table  1.  Plotting  Capability  Guide  (con't) 


StKnal  ratio-tXI 


(>ricfitation 


Service  voImc 


Signal  ratio  contours 


Sinilar  to  above,  but  abscissa  is  desired  racility-to- 
receiver  distance  and  the  station  separation  is  fixed. 

Undesired  station  antenna  orientation  witli  rcsitect  to 
the  desired  to  undcsired  station  line  versus  required 
facility  separation  curves  arc  plotted  for  several  de- 
sired station  antenna  orientations.  These  curves  shov 
the  naxiaua  se]’>aration  required  to  obtain  a specified 
D/U  signal  ratio  value  at  several  aircraft  locations 
(i.e.,  protection  points). 

Fixed  DA)  contours  arc  shown  in  the  altitude  versus 
distance  plane  for  a fixed  station  separation  and  tine 
availabilities  of  5,  50,  and  95  percent. 

Contours  for  several  D/U  values  are  sliown  in  the  alti- 
tude versus  distance  plane  for  a fixed  station  separa- 
tion and  tine  availability. 


* .Viditional  discussion,  by  capability,  is  provided  in  APTWCATICNS  OlIDE  [21,  sec^.  3.2).  _ 

»»  Applicable  only  to  the  line-of-slght  region  for  spherical  earth  geonetry.  Variability  with  tiae  and 
horizon  effects  are  neglected  and  the  counteipoise  option  is  not  available.  The  phase  change  a$so> 
dated  with  surface  reflection  in  the  lobing  region  is  token  as  0 or  1*0*  to  avoid  missing  Icfce  hulls. 


I 

1 

. taken  from  an  APPLICATIONS  GUIDIi  [21,  tables  1,  2,  3,  and  4]  where 

the  capabilities  and  input  requirements  are, discussed  in  detail. 
The  figure  numbers  in  table  1 refer  to  sample  capability  graphs 
contained  in  the  APPLICATIONS  CUTDn.  llcnce,  the  APPLICATIONS 
GlHDli  contains  one  or  more  sample  graphs  per  capability.  Even 
though  an  idea  of  the  capabilities  available  and  the  input  re- 
I quirements  can  be  obtained  from  these  tables,  anyone  seriously 

j considering  using  the  capabilities  -should  obtain  and  read  a copy 

! of  the  APPLICATIONS  GUIDE 

I This  repoi  t covers  exter.'' i one  that  were  made  to  IF-73  in  the 

process  of  develoring  the  1977  capabilities  of  table  1.  These 
extensions  allow  the  prog-  ' !'  to  be  used  for  a wider  variety  of 
problems  such  hs  rbose  ir..  elving  air/air  or  air/satellite  propa- 
j gation.  A brio  description  of  the  propagation  model  provided  in 
j section  2 is  followed  by  detailed  discussions  of  specific  model 
, extensions.  vii’i'O'-  change  * mid.'  iF-73  and  errata  for  the  1973 
; report  [17]  .ivo  c*''vercd  in  A]  p-.  rJix  A. 

Except  where  otherwise  indicated,  all  equations  provided 
here  are  dimensionally  consistent^  e.g.,  all  lengths  in  a partic- 
ular equation  are  expressed  in  the  same  units.  Calculations  are 
made  in  the  computer  programs  with  all  lengths  expressed  in  kilo- 
meters. Braces  are  used  around  parameter  dimensions  when  partic- 
ular units  are  called  for  or  when  a potential  dimension  diffi- 
culty exists.  A list  of  symbols  is  provided  in  Appendix  B. 


2.  PROPAGATION  MODEL 

The  IF-77  propagation  model  is  applicable  to  air/ground, 
air/air,  ground/satellite,  and  air/satellite  paths.  It  can  also 
be  used  for  ground/ground  paths  that  are  line-of-sight  or  smooth 
earth.  Model  applications  are  restricted  to  telecommunication 
links  operating  at  radio  frequencies  from  about  0.1  to  20  GHz 
with  antenna  heights  greater  than  1.5  ft  (0.5  m).  In  addition, 
radio-horizon  elevations  must  be  less  than  the  elevation  of  the 
higher  antenna.  The  radio  horizon  for  the  -higher  antenna  is 


Table  2.  Parameter  Specification,  General  [21,  table  2] 


^i!l!i!Liii^i'i?i-^i-'ii^ii.'iA'^'-J'ii<'jJ'i»''iiiiVji»'iCvj:i-i'»\~^^ 


taken  either  as  a common  horizon  with  the  lower  antenna  or  as  a 
smooth  earth  horizon  with  the  same  elevation  as  the  lower  antenna 
effective  reflecting  plane  [17,  sec.  A. 4.1;  21,  sec.  4.1],  Ran- 
ges for  other  parameters  associated  with  the  model  are  given  in 
table  2. 

At  0.1  to  20  GHz,  propagation  of  radio  energy  is  affected 
by  the  lower,  nonionized  atmosphere  (troposphere),  specifically 
by  variations  in  the  refractive  index  of  the  atmosphere  [2,  3,  4, 
8,  12,  20,  23,  35,  36].  Atmospheric  absorption  and  attenuation 
or  scattering  due  to  rain  become  important  at  SHF  [17,  sec.  A. 4. 5; 
23,  ch.  7;  36,  ch.  3;  42].  The  terrain  along  and  in  the  vicinity 
of  the  great  circle  path  between  transmitter  and  receiver  also 
plays  an  important  part.  In  this  frequency  range,  time  and  space 
variations  of  received  signal  and  interference  ratios  lend  them- 
selves readily  to  statistical  description  [20;  28;  31;  36,  sec. 
10].. 

Conceptually,  the  model  is  very  similar  to  the  Longley-Rice 
[26]  propagation  model  for  propagation  over  irregular  terrain, 
particularly  in  that  attenuation  versus  distance  curves  calcula- 
ted for  the  (a)  line-of-sight  [17,  sec.  A. 4. 2],  (b)  diffraction 
[17,  sec.  A. 4. 3],  and  (c)  scatter  (sec.  5)  regions  are  blended 
together  to  obtain  values  in  transition  regions.  In  addition,  the 
Longley-Rice  relationships  involving  the  terrain  parameter  Ah 
are  used  to  estimate  radio-horizon  parameters  when  such  informa- 
tion is  not  available  from  facility  siting  data  [17,  sec.  A. 4.1], 
The  model  includes  allowance  for 

(a)  average  ray  bending  [4,  (3.44),  (3.43),  (4.30);  5;  17, 
p.  44;  36,  sec.  4;  44] 

(b)  horizon  effects  [17,  sec.  A. 4,1], 

(c)  long-term  fading  [17,  sec.  A. 4;  36,  sec.  10], 

(d)  facility  antenna  patterns  [17,  sec.  A. 4. 2;  21,  sec. 

4.1], 

^The  numbers  in  parentheses  are  equation  numbers  for  the  given 
reference;  e.g.  [4]. 


i 


i 

I 


I 

I 

r 


I 


1: 


%' 


R: 


(e) 

(f) 
Cg) 

(h) 

(i) 
The 

[17,  sec 
(a) 

Cb) 

Cc) 


(d) 

(e) 

(f) 

(g) 

(h) 

(i) 

(j) 

Ck) 

(1) 

(m) 


surface  reflection  multipath  [6;  7;  16,  p.  17;  17,  sec. 
A. 6;  18,  sec.  CI-D.7], 

tropospheric  multipath  [3;  12,  sec.  3.1;  17,  sec.  A. 7; 
20;  25,  pp.  60,  B-2,  119], 

atmospheric  absorption  [15,  sec.  A, 3;  17,  sec.  A. 4. 5; 

36,  fig.  3.1], 

ionospheric  scintillations  [1;  16,  sec.  2.5;  18,  sec. 

evil;  32;  ,47],  and 

rain  attenuation  [11,  27]. 

IF-77  model  is  an  extended  version  of  the  IF-73  model 
A].  These  extensions  include  provisions  for 
sea  state  (discussion  of  this  extension  follows,  in  sec. 

3.1) , 

a divergence  factor  (sec.  3.2), 

a ray  length  factor  for  situations  where  the  free-space 
loss  associated  with  a surface  reflected  ray  may  be  sig- 
nificantly greater  than. that  associated  with  the  direct 
ray  (sec.  3,3), 

an  antenna  pattern  at  each  terminal  (sec.  3.4),  • 

circular  polarization  (sec.  3.5), 

frequency  and  temperature  variations  Of  the  complex  di- 
electric constant  for  water  (sec.  3.5),< 
long-term  power  fading  as  a function  of  time  block  (sec. 

4.2)  or  radio  climatic  region  (sec.  4.3), 
rain  attenuation  (sec.  4.4),, 
ionospheric  scintillation  (sec.  4.5), 

an  improved  method  for  calculating  the  transmission 
loss  associated  with  tropospheric  scatter  (sec.  5),; 
an  improved  estimate  of  the  distance  where  horizon  ef- 
fects can  be  neglected  (sec.  7), 

a free-space  loss  formulation  that  is  applicable  to  very 
high  antennas  (sec.  8), 

a formulation  for  facility  horizon  determinations  that 
includes  ray  tracing  (9.2), 
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(n)  ray  elevation  angle  adjustment  factors  to  allow  for  ray 

tracing  (sec.  10.2),  , 

(o)  antenna  tracking  options  (sec.  10.3),  and 

(p)  additional  antenna  pattern  options  [21,  pp.  85-88]. 

3.  EFFECTIVE  REFLECTION  COEFFICIENT 


The  formulations  used  previously  [17,  pp.  52-57  and  77-79] 
for  effective  reflection  coefficient  were  extended  to  permit 

(a)  surface  roughness  to  be  specified  by  sea  state  (sec.  3.1), 

(b)  both  antennas  to  be  high  (e.g.,  both  aircraft)  by  incor- 
porating allowances  for  a divergence  factor  (sec.  3.2) 
and  a ray  length  factor  (sec.  3.3), 

(c)  both  terminals  to  have  a vertical  antenna  pattern  asso- 
ciated with  them  by  using  a gain  factor  (sec.  3.4),  and 

(d)  circular  polarization  (sec.  3.5)  as  ah  option. 


3.1  Sea  State 

The  1977  computer  programs  allow  water  surface  roughness  to 
be  specified  by  sea  state  or  the  root-mean- square  (rras)  deviation, 

0.  , of  surface  excursions  within  the  limits  of  the  first  Fresnel 
h 

zone  in  the  dominant  reflecting  plane  [17,  p.  53;  26,  p.  3-23; 

36,  sec.  5.2.2].  Table  5 provides  the  relationship  between  sea 
state  and  that  is  used  in  the  model. 

Values  for  provided  in  table  5 were  estimated  using  sig- 
nificant wave  height,  estimates  from  Sheets  and  Boatwright 

[41,  table,!]  with  a formulation  given  by  Moskowitz  [29,  (1)]; 

1 . e . , 

where  and  have  the  same  units. 

Once  obtained,  values  of  are  used  as  they  were  in  IF-.73, 
to  calculate  "reflection  reduction  factors"  F^j^  r?^7»  (66,)^],,  and 
^ dah  C194)  as  corrected  in  Appendix  A of  this  repp,rt],.  iCpm- 

parisons  of  these  reflection  reduction  factor  formulations  with 
other  formulations  and  data  have  been  made  [18,  sec.  CI-D.7]. 


I 


I 


I 


I 

i 


I 
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3.2  Divergence  Factor 

The  divergence  factor,  D,  is  used  to  allow  for  the  divergence 
of  energy  reflected  from  a curved  surface  in  the  effective  reflec- 
tion coefficient  formulation.  It  is  defined  by  Reed  and  Russell 
[35,  p.  103]  "as  the  ratio  of  the  field  strength  obtained  after 
reflection  from  a spherical  surface  to  that  obtained  after  reflec- 
ting from  a plane  surface,  the  radiated  power*  total  axial  dis- 
tance, and  type  of  surface  being  the  same  i.n  both  cases,  and  the 
solid  angle  being  a small  elemental  angle  approaching  zero  in  mag- 
nitude." 

Figure  1 illustrates  the  geometry  for  reflection  from  a plane 
earth  and  a spherical  earth  where  the  relative  location  of  the 
source  reflecting  point  and  reference  plane  are  identical.  It 
also  shows  the  relative  size  of  the  ray  bundle  on  the  reference 
plane  for  each  case  (see  fig.  1.  caption).  The  divergence  factor 
is  related  to  the  reference  plane  area  associated  with  the  spher- 
ical earth  reflection,  and  the  plane  earth  reflection, 

by  

C2) 


D 


Derivations  of  expressions  for  D are  beyond  the  scope  of  this 
text,  but  such  developments  are  available  [7,  sec.  11.3;  8,  pp. 
DS-D?;  23,  sec.  5.2;  35,  sec.  4.27;  36],  An  exact  expression  for 
D that  is  very  similar  to  the  formula  provided  by  Beckmann  and 
Spizzichino  [7,  p.  223]  may  be  developed  by  extending  the  Riblet 
and  Barker  formulation  [37,  (13)]  to  the  special  case  where  prin- 
cipal radii  of  curvature  of  the  reflecting  surface  at  the  reflec- 
tion point  are  within,  a^,  and  normal,  a^,  to  the  plane  of  inci- 
dence. This  expression  is 


^2 

-1/2: 

2r,r~sinij/ 

(Tj  + rj)  Slni|» 

-1/2 


(3) 


'4 


16 


» 


p 


i 


where  the  ray  lengths  r^^  and  r^  along  with  the  grazing  angle  t|;  are 
shown  in  figure  2.  Here  the  r^^-  2**  and  a^  ^ must  be  expressed  in 
the  same  units;  e.g.,  kilometers.  For  the  spherical  earth  case, 


a^»a^*aj,,  so  that  (3)  may  be  expressed  as 


a 


D = 


1 + 


2R^(1  + sin'^  i<<) 


a sin  ij) 

a 


(^) 


-1/2 


f4) 


where 


C5) 


Values  for  ^ are  obtained  as  in  IF-73  [17,  p.  58],  and  values  for 
r^^  2 calculated  using 


1,2 


Hj^2  if  I*  - 90° 

otherwise 


(6) 


where  Hj^^2  defined  by  figure  2.  A formula  for  2 

is  included  Li  IF-73  [17,  p.  51].  * 

Divergence,  as  calculated  using  (4),  is  used  in  IF- 77,  It 
is  incorporated  into  IF-73  as  a factor  multiplying  the  left-hand 
side  of  equation  [17,  (68)]. 


3.3  Ray  Length  Factor 


The  ray  length  factor,  F^,  is  used  to  allow  for  situations 


where  the  free-space  path  loss  associated  with  the  reflected  ray 
may  be  significantly  gr,'ater  than  that,  associated  with  the  direct 
ray.  It  is  determined  using 

r 

(7) 


Ft. 

r r 
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'This  notation,  r^^  2>  Is  used  to  imply  r^  or  r2i 

18 


T?. 


where  is  the  direct  ray  length  and  rj^2  reflected  ray 

length  (rj^  + r2)  as  illustrated  in  figure  2.  Incorporation  into 
IF- 75  is  accomplished  by  using  it  as  a multiplying  factor  to  the 
left-hand  side  of  equation  [17,  (68)]. 


3.4  Gain  Factors 

The  antenna  gain  factors  g^  ^ and  are  used  to  allow 

for  situations  where  the  antenna  gains  effective  for  the  direct 
ray  path  differ  from  those  for  the  reflected  ray  path.  Figure  3 
illustrates  the  two-ray  path  and  indicates  the  gains  involved. 
These  are  the  relative  voltage  antenna  gains  (volts/volt  or  V/V) 


associated  with  the  direct  ray  at  terminal  one  or  two,  g^^^  2» 


those  associated  with  the  reflected  ray,  2*  They  are  measured 
relative  to  the  main  beam  of  their  respective  terminal  antenna; 
i.e.,  for  main  beam  conditions  gj^j^  2 " ^ This  conven- 

tion is  consistent  with  usage  in  IF-73  [17,  p.  39].  However,  it 
is  NOT  CONSISTENT  with  usage  in  the  Multipath  Handbook  [18,  sec. 
CI-D.3]  where  identical  symbols  are  used,  but  the  gains  are  mea- 
sured relative  to  an  isotropic  antenna. 

In  general,  these  gains  are  complex  quantities,  but  IF- 77  in- 
cludes provisions  for  scalar  gains  only.  In  many  practical  appli- 


elevation  angles  .where  the  relative  phase  is  either  expected  to 
be  near  zero  or  is  unknown,  so  that  the  complex  nature  p£  these^ 
gains  is  largely  academic.  They  are  called  voltage  gains  .since 
they  are  a voltage  ratio  that  could  be  considered  dimensipnless 
(volt/volt);,  but  are  different  from  gains  expressed  as  power  ra- 
tios (watt/watt)  that  could  also  be  considered  dimensionless.  De- 
cibel gains  above,  main  beam  values  are  related  to  these  gains  by  - 
formulas  such  as  • 

(8)‘ 


and 


|gRl_2ltV/Vl 

The  formulations  for  gjj  are 

1|8d1  ®D2  linear  polarization 


2d[V/V3 


and 


8r[V/V1 


<=hD2^%Dl  *VD2> 


/I  for  omnidirectional  antennas 
and/or  circular  polarization 
( see  text  below) 


8r1  ^R2 


otherwise 


(9) 


(10) 


(11) 


where  omnidirectional  implies  that,  for  the  radiation  .angles,.  6£ 
interest,  gj^j^  ■ ®nd  gj^2  “ 8i)2*  problems  involving,  circular 
polarization,  horizontally  polarized  (g^jjj^  2 ^hRl  2^ 
vertically  ipolarize.d  (g,yjjj^  2 8vri  2^  99nippn.ejits  are  j^^Ad. 
Linear  polarization  is  considered,  to  be  either  vertical  Ip.r  hori- 
zontal with  the  polarization  associated  wLth^.  .gp.. selected  ; , j 

accordingly.  Defining  g^  as  1 for  circular  polarization  is  done 
to  allow  the  antenna  gains  to  be  included  in  the.  reflect  ion 
coefficient  formulation  of  lE-T/  in  a simple. Kay  for  horizontal  or 
vertical  polarization.  The  circular  polarization  case 

will  be  discussed  in  the  next  section. 
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The  IF- 73  was  extended  to  include  gj^  by  incorporating  it  as 
a replacement  for  the  multiplying  factor  g in  left-hand  side  ofr 
two  equations  [17,  (68) j (69)].  For  special  cases  where  only  the 
antenna  patterns  of  IF-73  are  used  (i.e.,  isotropic  aircraft  an- 
tenna pattern  and  specific  facility  antenna  patterns) , gj^  reduces 
to  g [17,  (67)]  . The  effects  of  g^  are  included  in  IF-73  [1-7, 
(8i),  (82)]  with  a variable  named  g^,  but  since  gjj  can  now  be 
complex,  it  is  necessary  to  use  |gjjl  in  one  of  the  IF-73  equations 
[17,  (82)].  In  addition,  it  should  be  realized  that  the  aircraft 
antenna  gain  is  not  necessarily  0 dBi  as  it  was  in  IF-73  [17, 
p.  37]. 

The  gain  factor  gj^^  is  similar  to  gj^  except  that  gj^^  involves 
gains  8yj^2  2’  i«e.. 


Also 


gRyfV/V] 


1 1 for  omnidirectional  antennas 


®vRl  ®vR2 


ptherwise 


% 


(12) 


8Rh(V/Vl 


1 for  omnidirectional  antennas 


(13) 


( ®hRl  ®hR2  otherwise  | 

where  gj^j^  is  for  horizontal  polarization.  These  factors  will  be 
used  in  the  formulation  of  complex  plane  earth  reflection  co- 
efficients for  circular  polarization  that  is  given  in  the  hexf 
section. 

3.5  Plane  Earth  Reflection  Coefficient 
Values  for  the  "complex  plane  earth  reflection  coefficient , - 
R exp(-j*),  used  in  IF-73  [17,  pp.  52,  53]  depend  ‘on  the  rela- 
tive dielectric  constant,  e,  and  conductivity,  a,  along' with  ' 
wavelength,  x,  grazing  angle,  ip,  and  polarization  [7,  p.'  219; 

18,  sec.  CI-D.8;  23,  p.  396;  35,  p.  88;  36;  sec’.  Ili.lj.  For' 
Vertical  polarization  (electric  field  in  thk  plane  of  iricfiderici^i 


or  horizontal  polarization  (electric  field  normal  to  plane  of 
incidence)  R exp(-j4i)  is  given  by  


Ry  ®^p['  j “ ®v  J 


e^sin(if/)-Y 


vr  ®R 


respectively,  where 


j-  1 sin(ij<)-Y 

*h  "Pt  j ’■"  - '=h>J  - sin(i|.)*Y,  ®R’ 


Y =We  - cos  t|) 
c f c 


is  complex,  the  complex  relative  dielectric  constant,  e^,  is  de- 
fined as 

•e^  * e-j  60X  0 , .,(17) 

and  gj^  is  from  (11) . 

In  IF-77,  linear  polarization  gain  factors  (sec.  3.4)  and  re-' 
flection  coefficients  are  combined  to  obtain  a reflection  coeffi- 
cient formulation  for  circular  polarization;  i.e.. 


'«Rv*v  • 


'«Rv*v  • 

This  formulation  is  used  only  for  antennas  with  the  same  polari- 
zation sense  (e.g.,  both  right-handed). 

For  a perfect  dielectric  (a  * 0 so  that  e^*e),  the  numerator 
of  (14)  win  go  to  zero  when  'I' * ipg  where. 


<(n  * Sin 


._-l 


so  that  Ry-0.  This  critical  angle  is  called  the  Brewster  angle 
and  a silnilar  angle  associated  with  reflection  from  a surface 
that  has  non- zero  conductivity  is  called  the  pseudo^Brewster  angle 
[36,  sec.  III.l].  Equation  (19)<  may  be  used  to  estimate  the 
pseudo  Brewster  angle  when  e > 60Aa.  Figure  4 shows  the  dip 


Figure  4.  Comparison  of  reflection  coefficients  for  sea  water 
vertical  polarization  [18,  p.  CI-103], 


in  the  reflection  coefficient  for  vertical  polarization  associa- 
ted with  the  pseudo-Brewster  angle  along  with  the  abrupt  change 
in  phase  that  occurs  as  ip  goes  through  its  critical  value.  This 
change  in  phase,  which  does  not  occur  for  horizontal  polarization, 
will  change  the  rotation  sense  of  circularly  polarized  waves  that 
are  reflected  from  the  surface;  i.e.,  when  a circularly  polarized 
wave  is  reflected,  its  rotation  sense  will  remain  unchanged  only 
if  the  grazing  angle  is  less  than  the  pseudo-Brewster  angle. 

In  IF- 77,  e and  a for  water  may  be  estimated  with 


s s o 


1 + C2irfT) 


+ e 
2 o 


o [mho/m]  « f^T(e  - ejj)/2865  + (21) 

where  e„  is  the  static  dielectric  constant,  = 4.9  is  the  di- 
s o 

electric  constant  representing  the  sum  of  electronic  and  atomic 
polarizations,  f[MHz]  is  frequency,  T[ps]  is  relaxation  time,  and 
(mho/m]  is  the  ionic  conductivity.  Values  of  e^,  T,  and  ob- 
tained using  Saxton  and  Lane  [40]  are  provided  in  table  6 for 
fresh  water  and  sea  water. 

Figure  4 provides  a comparison  of  reflection  coefficients  cal- 
culated using  the  fixed  values  (table  6)  of  surface  constant  given 
by  Rice  et  al.  [36,  p.  III-7]  for  sea  water  (solid  lines)  with 
those  determined  via  calculated  surface  constants  (dashed  lines). 
More  such  comparisons  are  available  [18,  sec.  CI-D.8]. 

4.  VARIABILITY 

Model  extensions  that  are  concerned  directly  with  transmis- 
sion loss  (or  received  signal  level)  variability  are  discussed 
this  section.  These  extensions  Include  pro visions  for  (1)  mixing, 
distributions  (sec.  4.1),  (2)  computing  long-term  variability'  for 
various  time  blocks  (sec.  4.2)  or  climates  (sec.  4.3>,  aihtf  i3)' 
estimating  the  effects  of’ rain  attenuation  (sec.  4.4):  and'  loho- 
spheric  scintiHatibn-  (sec.  -AVSi.  ' ' 
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Table  6.  Surface  Types  and  Nominal  Constants 


Surface  type 

e 

0 [mho/m] 

Poor  Ground 

4 

0.001 

fa") 

Average  Ground  '■  •' 

15 

0.005 

Good  Ground 

25 

0.02 

Fresh  Water 

81 

0.01 

Sea  Water 

81 

5 

rbi 

Concrete  *•  ^ 

5 

0.01 

Metal 

10 

107 

For  Fresh  Water 

0°C 

10“C  20“C 

88 

84  80 

1.87 

X 10"5 

1.36  X 10-5  i.oi 

[mho/m] 

0.01 

0.01  0.01 

s 

For 

Sea  Water 

O^C 

lO^C  20‘’C 

75 

72  69 

T[ps] 

1.69 

X 10-5 

1.21  X 10-5  9.2  : 

0^  [mho/m] 

3.0 

4.1  5.2 

^^^From  Longley  and  Rice  [26, 
Estimated  [2'i,  p.  398]. 
Estimated  [34,  p.  235,  p. 
^‘^^From  Saxton  and  Lane  for  0 
^®^From  Saxton  and  Lane  for  3 

table  2]. 

240] . 

5 salinity  [4.0,  table  1]. 

.6%  salinity  [40,-  table  1] . . 
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4.1  Mixing  Distributions 

Subroutines  have  been  incorporated  into  the  computer  programs 
to  allow  the  distributions  that  characterize  portions  of  the  vafif 
ability  associated  with  a particular  model  component  to  be  mixe'd' 
in  order  to  obtain  the  total  variability  for  that  component.  For^ 
example,  different  fractions  of  the  time  may  be  characterized  by 
signal  level  distributions  associated  with  different  Ionospheric 
scintillation  groups,  and,  with  these  subroutines,  they  can  be 
weighted  and  combined  (mixed)  to  obtain  the  total  variability  as- 
sociated with  ionospheric  scintillations  (sec.  4.5). 

The  process  of  mixing  N cumulative  variability  distributions 
may  be  summarized  as  follows: 

1)  Select  M (ten  or  more)  levels  of  variability 

Vj^,  ...,  ...»  that  cover  the  entire 

.range  of  the  transmission  loss (or  power  available,  etc.) 
values  involved. 

2)  Determine  the  fraction  of  time  (weighting 
factor)  for  which  each  distribution  is 


W. 


and  - 


applicable;  i.e.,  W^,  . . . , , ...,  ,.j^. 

3)  Determine  the  time  availability  (fraction  of 

time  during  which  a distribution  is  appli- 
cable that  a specific  level  of  transmission 
loss  is  not  exceeded)  for  each  distribution 
at  the  selected  levels;  i.e.,  ..., 

‘^ij’  ^MN- 

4)  Calculate  time  availabilities  for  the  mixed 
distribution  that  corresponds  to  the  variability 
levels  selected,  i.e.. 


‘^l  ' ^11  ''^l  * 


‘li^^il^l^  •••  ■^‘liN^ 


‘Im  " ‘^Ml  "^1  • • • + ‘?Mj  <1mN 
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This  process  is  the  same  as  the  one  used  by  Rice  et  al. 

[36,  sec.  Ili.7.2]  to  combine  transmission  loss  distributions 
for  time  blocks  (sec.  4.2)  to  obtain  distributions  for  summer 
and  winter.  It  is  also  essentially  the  same  as  the  method 
recommended  by  Whitney  et  al.  [46,  p.  1099;  47,  sec.  6]  to  com- 
bine distributions  of  fading  associated  with  various  ionospheric 
scintillation  index  groups  (sec.  4.5). 

When  this  process  is  used  to  mix  distributions  of  long-term 
variability,  the  required  variability  functions  are  obtained  from 

V^.(q)  »V(0.S)  +Y(q)|^  (23) 

where  indicates  that  the  V(0.5)  and  Y(q)  are  appropriate  for 

the  conditions  (time  block  or  climate)  associated  with  a particu- 
lar value  of  the  subscript  c.  For  example,  V(0.5)  and  Y(q)  values 
for  different  climates  can  be  obtained  with  the  information  sup- 
plied in  section  4.3,  and  mixing  can  be  used  to  estimate  variability 
for  areas  near  a border  between  two  different  climate  types.  Af- 
ter mixing,  Y(q)  values  needed  for  later  calculations  may  be  ob- 
tained from  using 

• Y(q)  »V(q)  - V(0.5)  (24) 

where  all  variables  in  (24)  are  associated  with  the  resulting 
mixed  distribution.  Similarly,  when  mixing  variabilities  associ- 
ated with  ionospheric  scintillation, 

Yi,(q)-YjC,)|c,  (25, 

and  the  distribution  resulting  from  the  mixing  is  taken  as  Yj(q) 
for  later  calculations. 

4.2  Time  Blocks 

The  long-term  variability  portion  of  IF-73  [17,  sec.  4.5] 
has  been  extended  to  allow  the  variability  associated  with  speci- 
fic time  blocks  or  a combination  of  time  blocks  (sec.  4.1)  to  be 
used.  Table  7 shows,  the  months  and  hours  of  the  day  that  cor- 
respond to  the  various  time  blocks,  these  blocks  and  season 
groupings  are  used  to  describe  the  diurnal  and  seasonal  varia- 
bility in  a continental  temperate  climate  [36,  sec.  III. 7.1]. 
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Wble  7. 

Time  Block  Ranges  [36,  sec. 

ni.7.1]. 

No. 

Months 

Hours 

1 

Nov. 

- Apr, 

0600  - 1300 

2 

Nov. 

- Apr. 

1300  - 1800 

3 

Nov. 

- i^r. 

1800  - 2400 

4 

May 

- Get. 

0600  - 1300 

5 

May 

- Oct. 

1300.  - 1800 

6 

May 

- Oct. 

1800  - 2400 

7 

May 

- Oct. 

0000  - 0600 

8 

Nov. 

- .^r. 

0000  - 0600 

Sumner 

May 

- Oct. 

all-hours 

Winter 

Nov. 

- .^r. 

all-hours 

Variability  associated  with  the  time  blocks  and  seasons  given  in 
table  4 were  incorporated  into  IF- 77  by  allowing  the  constants 
given  by  Rice  et  al.  [36,  tables  III, 2,  III, 3,  and  III. 4]  to  be 
used  in  an  equation  of  IF- 73  [17,  (178)], 

4.3  Climates 

The  IF-77  includes  extensions  that  allow  the  use  ,of  long-term 
power  fading  (variability)  -applicable  to  various  climates..  How- 
ever, the  long-term  variability  of  IF-73  [17,  sec.  A. 5]  is 
normally  used  in  IF-77  except  when  another  climate,  time  block;,, 
(sec.  4.2)  or  combination  (sec.  4,1)  of  climates  (or  time  blocks) 
is  specifically  requested.  The  ability  to  mix  distributions 
(sec.  4.1)  that  characterize. long-term,  power  fading  to  obtain 
combinations  of  climates  or  time  blocks  [36,  sec.  III.7.2]  adds  i 
flexibility  to  the  model.  For  example,  (a)  more  meaningful 
comparisons  can  be  made  with  data  in  cases  where  the  data  col- 
lected do  not  represent  all  hours  of  the  day  or  all  months,  of 
the  year  [14,  sec.  4.3],  and  (b)  variability  formulations 

that  may  become  available  for  prop,agation  via  specific  mechanisms 
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(forward  scatter,  diffraction,  partial  reflections,  ducting,  etc.) 
can  be  comliined  in  accordance  with  the  fraction  of  the  total  time 
that  they  are  effective. 

The  various  climate  types  are  listed  in  table  8 , including 
supplementary  data  to  aid  in  the  selection  of  the  appropriate 
type  for  a specific  radio  link.  Table  8 (Samson  and  Hart,  DOC-BL, 
informal  communication)  is  based  primarily  on  the  annex  to  CCIR 
Report  244-2  [10],  and  is  presented  here  as  the  best  available 
information  in  lieu  of  maps.  If  a path  is  near  a border  between 
two  different  climate  types,  calculations  can  be  made  for  each 
climate  or  mixing  (sec.  4.1)  can  be  performed  to  combine  the 
variabilities  associated  with  the  climates  involved. 

The  formulation  for  long-term  variability  given  here  as  a 
function  of  effective  distance,  d^  [17,  (177)  on  p.  751,  is  based 
on  curves  provided  in  CCIR  Report  244-2  [10].  Algebraic  expres- 
sions fitted  to  the  modified  versions  of  the  CCIR  curves  are  used 
(Hufford  and  Longley,  DOC-BL,  informal  communication).  It  was  felt 
that  the  CCIR  estimates  for  Climates  3,  7a,  and  7b  are  not  typical 
for  longer  distances  and  values  of  time  availability  of  1 percent 
or  less;  i.e.,  time  fraction  of  q_<0.01.  The  near  free-space  values 
shown  by  the  CCIR  curves  for  paths  with  d^>400  km  require  that 
the  signal  be  carried  within  a duct,  and  while  this  could  occur, 
it  is  hot  considered  typical  enough  to  be  included  in  a general 
variability  formulation.  Thus,  curves  developed  from  the  formu- 
lation provided  here  would  differ  somewhat  from  the  applicable 
CCIR  recommendations  and  reports,  but  they  are  thought  to  be  im- 
proved estimates.  A more  complete  discussion  of  this  formulation 
that  includes  graphs  for  various  climates  has  been  prepared  for 
publication  in  a Military  Handbook  (MIL  HDBK  417)  titled  ’’Facility 
Design  Handbook  for  Transhorizon  Communications”. 

The  formulation  is  incorporated  into  computer  programs  via 
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V(0.5) 

•Yo(0.9) 


{ 

l+(dg/bi)2 

lM(d^-b2)/b3l2 

« 

55^ 


(26) 


and 


Y(O.l)  - Yo(O.l)  g(0.1,f)  (27) 

YC0.9)  - Yq(0.9)  g(0.9,f)  (28) 

where  values  obtained  are  used  as  in  IF- 73  for  V(0.5)  [17,  (190)], 
Y(0.1)[17,  (180)],  and  Y(0.9)  [17,  (181)].  liffective  distance. 


d , is  determined  as  it  was  in  IF-73  [17,  (177)];  values  for  the 
6 


constants  b^,  b2,  bj. 


Cj^,  and  C2  to  be  used  for  each  climate  are 


provided  in  table  9 and  the  factors  g(0.1,£)  and  g(0.9,f)  are 
calculated  as  follows: 


g(0. !,£)--< 


1 for  all  climates  except  2,4,  and  6, 

0.18  sin  5 logjQ(f/200)+  1.06 

for  60  <.fi  1500  MHz  in  Climates  2 and  6, 
1 suggested  for  60  <£<  200  MHz  in  Climate  4, 
0.10  sin  5 logj^Q(f/200)+  1.02, 

for  200  _<f  < 1500  MHz  in  Climate  4, 

0.93  for  f >1500  MHz  in  Climates  2,4,  and  6 


y (29) 


and 


g(0.9,f) 


1 for  all  climates  except  6, 

0.13  sin[S  logjQ(f/200)]+  1.04 

tor  60££i  1500  MHz  in  Climate  6, 
0.92  for  f > 1500  MHz  in  Climate  6 


(30) 


Note  that  the  above  formulation  is  incomplete  in  some  respects, 
but  that  approximations  are  suggested  to  fill  the  gaps;  i.e., 
(a)  g(0.1,f)  in  (29)  is  approximated  by  1 for  60  ^£<200  MHz 
in  Climate  4,  and  (b)  the  constants  for  Climate  8 (table  9) 
are  approximated  with  those  of  Climate  6. 
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Table  9.  Constants  Used  to  Calculate  V(0,5),  Y^(0,1),  and  Y^  CO^S) 


Climate 

Parameter 

>>1 

•^2 

>>3 

^1 

^2 

1.  Equatorial 

V(0.5) 

144.9 

190.3 

133.8 

-9.67 

12.7 

636.9 

134.8 

95.6 

2.70 

131.1 

Yq(0.9) 

762.2 

123.6 

94.5 

-2.73 

-204.4 

2.  Ccsitinental 

V(0.5) 

228.9 

205.2 

143.6 

-0.62 

9.19 

subtropical 

y^co.!) 

138.7 

143.7 

98.6 

8.8 

19.9 

Y^co.g) 

100.4 

172.5 

136.4 

-3.41 

-9.83 

3.  Maritime 

V(O.S) 

262.6 

185.2 

99.8 

1.26 

15.5 

subtropical 

Y^(O.l) 

165. 3 

225.7 

129.7 

12.9 

12.3 

Yq(0.9) 

138.2 

242.2 

178.6 

-7.83 

-8.52 

4.  Desert 

V(0.5) 

84.1 

101.1 

98.6 

-9.21 

9.05 

Y^(oa) 

464.4 

93.1 

94.2 

4.72 

204.2 

Y^(0.9) 

139.1 

132.7 

193.5 

-2.54 

-16.8 

5.  ^Mediterranean 

6.  Continental 

V(0.5) 

228.9 

205.2 

143.6 

-0.62 

9.19 

ten5)erate 

Y^(O.l) 

93.2 

135.9 

113.4 

6.04 

10.4 

Y^(0.9) 

93.7 

186.8 

133.5 

-3.43 

-9.17 

7a.  Maritime 

VC0.5) 

141.7 

315.9 

167.4 

-0.39 

2.86 

temperate 

Y^(O.l) 

216.0 

152.0 

122.7 

11.0 

17.9 

overland 

Yq(0.9) 

187.8 

169.6 

108.9 

-8.79 

-13,3 

7b.  Maritime 

V(0.5) 

2222.0 

164.8 

116.3 

3.15 

857.9 

temperate 

Yo(O.l) 

136.2 

188.5 

122.9 

10.8 

10.5 

oversea 

Y^CO.9) 

609.8 

119.9 

106.6 

-10.9 

-217.6 

8,  *Polar 

Use  climate  6 

* For  climates  numbers  5 and *8,  Mediterranean  and  Polar, values  are  not 
available:  a substitute  for  Polar  is  suggested  for  use  unless  more 
definite  information  is  available  from  other  sources. 
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4.4  Rain  Attenuation 

The  rain  attenuation  model  used  in  IF-77  is  largely  based  on 
material  in  informal  papers  by  C.  A.  Samson  (DOC-BL)  on  "Radio 
propagation  through  precipitation"  and  "Rain  rate  distribution 
curves".  Only  those  portions  of  these  papers  that  are  directly 
related  to  IF-77  are  repeated  here.  However,  an  attempt  has  been 
made  to  cite  references  on  which  his  work  is  based. 

Two  options  for  rain  attenuation  are  available  in  IF-77.  The 
first  is  for  use  in  a "worse  case"  type  analysis  where  a particu- 
lar rainfall  attenuation  rate  is  assumed  for  the  in-storm  path 
length,  and  the  additional  path  attenuation  associated  with  rain 
is  simply  taken  as  the  product  of  this  attenuation  rate  (in  dB/km) 
and  the  in-storm  ray  length.  This  ray  length  is  determined  in  ac- 
cordance with  the  method  discussed  as  step  4 of  option  two. 

Option  two  involves  computer  inputs  of  rain  zone  (which  de- 
termines a rainfall  rate  distribution)  and  storm  size.  Rain  zones 
may  be  estimated  using  figure  5 or  6,  and  the  storm  size  (diameter 
or  long  dimension)  is  assumed  to  be  one  of  three  options:  5,  10, 
or  20  km  (corresponding  approximately  to  a relatively  small,  aver- 
age, or  very  large  thunderstorm).  The  maximum  distance  used  in 
calculating  path  attenuation  with  this  option  is  the  storm  size 
since  it  is  assumed  that  only  one  storm  is  on  the  path  at  a time. 
The  process  used  to  include  rain  attenuation  estimates  in  IF-77 
for  this  option  may  be  summarized  as  follows: 

1)  Determine  point  rain  rates . Point  rain  rates 
(rate  at  a particular  point  of  observation) 
not  exceeded  for  specific  fractions  of  the 
time  are  determined  from  table  12  for  the 
rain  zone  of  interest.  Values  listed  in  this 
table  were  taken  from  estimated  distributions 
[22;  38;  39;  45]. 
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Table  12. 

Point  Rain  Rates 

(mm/hr)  not  Exceeded  for  a 

Fraction  of  Time,  q. 

Rain  Zone 

r 

1 

2 

3 

4 

5 

6 

<.98 

0 

0 

0 

0 

0 

0 

.99 

0.17 

0.25 

0.31 

0.54 

0.75 

1.0 

.995 

0.62 

0.98 

1.54 

2.07 

2.7 

3.35 

.998 

1.8 

3.1 

4.8 

6.2 

7.8 

9.4 

.999 

3.2 

5.4 

8.8 

11.7 

14.0 

17.0 

.9995 

5.1 

9.6 

14.5 

19.0 

23.5 

28.5 

.9998 

8.2 

17.0 

25.0 

33.0 

40.0 

48.0 

.9999 

11.3 

22.8 

34.0 

44.5 

54.0 

67.0 

.99995 

14.6 

29.5 

43.0 

57.0 

68.0 

84.0 

.99998 

18.8 

37.8 

56.0 

73.0 

91.0 

112.0 

.99999 

24.0 

44.0 

64.0 

86.0 

110.0 

160.0 

2)  Determine  path  average  rain  rates.  Each  point  rain 
rate  resulting  from  step  1 is  converted  to  a path 
average  rain  rate  by  using  linear  interpolation  to 
obtain  a multiplying  factor  from  the  values  provided 
in  table  13.  These  values  were  taken  from  curves 
fitted  to  data  collected  in  Florida  [22]. 

3)  Determine  attenuation  rate.  For  each  path  average 
rain  rate  resulting  from  step  2,  an  attenuation  rate 
A^^(q)  [dB/km]  is  determined  using  linear  interpolation 
between  the  values  provided  in  table  14.  These  are 
theoretical  values  [27]  that  were  determined  for  a 
Laws  and  Parsons  [24]  drop  size  distribution. 

4)  Determine  the  in- storm  ray  length.  First  the 


length  of  the  direct  ray  that  is  within  T^^ 
of  the  earth's  surface  is  determined  using  the 
methods  described  in  1973  for  the  calculation 


^eo,w  replace 

•'“’'  ^o.5,»  ’’eo.w  Teo,s,w'-  ^es  ' 

is  taken  as  the  storm  size;  i.e.,  storm  height  is 

assumed  to  be  equivalent  to  storm  diameter  [9,  p. 

98].  Then  the  final  in-storm  ray  length,  r^,  is 

calculated  using 


r^ [km] » 


IT  if  r > T 
es  es-  esi 


fr,  . otherwise 
6 :> 


(33) 


5) 


For  transhorizon  paths,  the  storm  is  assumed  to  be 

between  the  facility  and  its  horizon  so  that  r is 

es 

not  increased  because  of  ray  lengths  within  T of  the 

G S 

surface  that  occur  beyond  the  facility  horizon. 
Determine  rain  attenuation  values . Values  for  the 
attenuation,  Aj.(q)  for  a particular  fraction  of  time 
are  calculated  using 


Aj.(q)[dB] 


0 for  q<.0.98 
A^y(q)rg  otherwise! 


(34) 


where  Ay^(q)  values  come  from  step  3 and  the  value 
for  r^  is  from  step  4.  Note  that  table  12  yields 
distributions  of  rain  attenuation  that  are  zero  for 
q < 0.98 . 

6)  Combine  rain  attenuation  variability  with  other 
variabilities . Variability  for  rain  attenuation 
Yj.(q)  is  related  to  the  distribution  of  rain  attenu- 
ation from  (34)  by 

Yr(q)  - -Aj.(q)  . (35) 

It  is  combined  with  the  long-term  power  fading  vari- 
ability, Yg(q)  [17,  sec.  A. 5],  and  multipath  variability, 
Yjj(q)  [17,  p.  38],  by  i^luding  Y^^Cq)  in  an  equation 
of  IF-73  [17,  (5)];  i.e., 
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Vjtq)  - i /v^2((,.),*Yj,2(q)  *Y^2(q)  +Yj^Cq)  dB  (36) 

+ for  q ^ 0, 5 
- otherwise 

where  total  variability  and  Yj(q)  is  a 

variability  included  in  IF- 77  to  allow  for  ionospheric 
scintillation  (sec.  4.5). 
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4.5  Ionospheric  Scintillation 

Variability  associated  with  ionospheric  scintillation  [1;  46] 
for  paths  that  pass  through  the  ionosphere  (i.e.,  on  earth/satel- 
lite paths)  at  an  altitude  of  about  350  km  [32,  p.  4]  is  included 
in  IF-77.  This  variability,  Yj (q)  dB,  is  determined  using  figure 
7 [46;  47,  fig.  6]  directly  if  calculations  are  to  be  for  a spe- 
cific scintillation  index  group  (see  fig.  7 inset)  or  using  a 
weighted  mixture  of  the  figure  7 distributions  (sec.  4.1)  where 
the  weighting  factors  are  estimated  for  specific  problems.  For 
example,  a computer  program  available  at  NTIA/ITS  [33]  that  is  an 
extension  of  the  Fremouw  model  [13]  can  be  used  to  estimate  weight- 
ing factors  for  frequencies  up  to  400  MHz  [32].  An  equation  given 
previously  in  section  4.4,  (36),  is  used  to  add  Yj(q)  to  IF-77. 

Provisions  exist  (table  2,  index  group  6)  to  allow  Yj(q) 
to  change  with  earth  facility  latitude  when  a geostationary 
satellite  is  involved  and  the  earth  facility  locations  are  along 
the  subsatellite  meridian.  Figure  8 shows  the  distributions  cur- 
rently used  when  this  option  is  selected.  These  distributions 
were  developed  by  mixing  distributions  for  particular  scintilla- 
tion index  groups  in  accordance  with  the  estimated  time  for  which 
they  would  b"*  present  at  a frequency  of  136  MHz  so  that  the  fre- 
quency scaling  factor  discussed  below  should  be  used  with  these 
distributions  [43,  table  5].  However,  only  minor  program  modi- 
fications would  be  necessary  to  incorporate  other  distributions 
that  might  be  of  interest. 
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When  the  distributions  of  figure  8 are  used  for  a frequency 
other  than  136  MHz,  an  optional  frequency  scaling  factor  should 
be  used.  It  relates  Y^Cq)  to  Y^^j^Cq)  from  figure  8 by 


Yj(q)  = (136 /f)"  Y^3g(q) 


(37) 


where  n varies  with  earth  facility  latitude,  epj^[43,  (27)];  i.e., 
1 for  0pj^  < 17“  or  9pj^  > 52“  , 

1 + (0pj^-  17)/7  for  17“<epL<24“ 


n = 


(38) 


2 for  24“<e^,<45“ 
FL 


1 + (52-0pL)/7  for  45“<0pL<52. 

Other  scaling  factors  could  be  used  with  minor  program  modifica- 
tions. 


5 . TROPO.SPHERIC  SCATTRR 

The  Rice  et  al.  [36,  sec.  9]  method,  which  is  used  to  calcu- 
late attenuation  for  tropospheric  scatter  in  IF-73  [17,  sec.  A. 4. 4], 
is  not  applicable  to  paths  that  involve  a very  high  antenna  such 
as  a satellite.  This  method  was  reformulated  by  Dr.  George  A.  Huf- 
ford  (DOC-BL,  informal  communication)  to  include  geometric  para- 
meters associated  with  very  high  antennas  where  these  parameters 
are  determined  using  ray  tracing  techniques.  The  resulting  formu- 
lation has  been  incorporated  into  IF- 77  and  is  presented  here. 

It  was  developed  using  kilometers  as  a measure  of  length  so  that 
all  lengths  in  the  formulas  of  this  section  are  in  kilometers. 

Frequency  and  the  basic  geometric  configuration  for  the 
tropospheric  scatter  path,  are  assumed  to  be  known  so  that  values 
for  the  following  parameters  are  available  where 

hi, 2 [km-msl]  * antenna  elevations  above  mean  sea 
level  (msl) , 

^rs  ® effective  reflecting  surface  elevation 

above  msl. 
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where 

h^g  is  also  taken  a:,  average  terrain  elevation, 

h [km-hrs]  = common  volume  elevation,  above  b^  , 

n,  2[bm]  = antenna  to  common  volume  ray 

Ng  [N-units]  = surface  refractivity  at 
0[rad]  = scattering  angle, 

lengths , 

^rs’ 

and 

X [km]  = wavelength. 

The  scattering  efficiency  term,  S , is  calculated  as  follows: 

t,  = 0. 031-2. 32(10'-^)N^  + 5.67(10‘°)N‘^, 
i s --  s 

(39) 

Y = 0.1424{l+c^exp[-(0.25h^)°]}, 

(40) 

e,  = 0.002Ng^-0.06Ng  + 6.6, 

(41) 

and 

p 

S=91.1-  ^ ^ 

® l+0,7716h^^ 

+ 20  log[(0.1424/Y)^exp(Yh  )]  . 

(42) 

The  scattering  volume  term,  S , is  calculated  as  follows: 

^1-^2 
s = ^ 

(43) 

*'l  + *^2 

where 

5 is  the  modules  of  asymmetry, 

A ri  2.2 

A = (l-s  ) , 

(44) 

i = 

(45) 

where 

1 is  the  total  ray  lengt). , 

n = yQl/2 , 

(46) 

Xj^  = (l+s)^n, 

(47) 

X2  = (l-s)^n. 

(48) 

K =*  2Tr/A 

(49) 

where 

K is  the  wave  number, 

^1,2  ' ^''®*^^l,2’^rs^  ’ 

(50) 

^^1,2  ""  ^^l,2'^’^^^‘^Pl,2’ 

(51) 
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B = 6 + 8s‘ 


(52) 


+ 8{l  + s)^2^P2^/q^, 

) ■)  ) O 

««il  '1  ■ >•  I'f  j'/'l  j’ 

+ 2(l-s)^l  + 2xj^^/qj^)(l+2x2/q2)» 


and 


C = 12* 


Sy  = 10  log 


/pi^/2\^  /P2*'^\^  P1^P2 

V P,  / V P2  ; Pi-P2^ 


2 2 

Pi  P2 


z/r 


(53) 


(54) 


Finally,  the  attenuation  for  scatter,  A , relative  to  free 


space  is  calculated  using 


AjdB]  =Sg  + Sy  + 10  log  0^  , 


(55) 


Xt 


6.  CONDITIONAL  ADJUSTMENT  FACTOR 


The  cond:. -'ional  adjustment  factor,  Ayj  is  used  in  IF- 73 

to  prevent  available  signal  powers  from  exceeding  levels 
expected  for  free-space  propagation  by  unrealistic  amounts  when 

the  variability  is  large  and  the  calculated  leference  level  is 
near  its  free-space  value.  This  is  accomplished  by  adding  Ay 
to  calculated  reference  basic  transmission  loss,  L^^.,  in  the 
computation  of  median  basic  transmission  loss,  (0.5)  [17,  pp. 
40,  41j.  However,  the  resulting  increase  in  transmission  loss 
can  be  too  large  for  frequencies  near  400  MHz  when  climates  or 
time  blocks  with  large  variabilities  are  used.  For  example, 
the  use  of  Time  Block  7 (sec.  4.2)  can  result  in  an  A^  of  20  dB. 

To  prevent  excessive  loss  increases  associated  with  Ay,  a 
formulation  to  keep  Ay^lO  dB  has  been  incorporated  into  IF- 77* 
This  formulation  may  be  summarized  as  follows; 
f 


■©h 

^bf 


= elevation  angle  correction  factor  [17,  (179)], 

» basic  transmission  loss  for  free  space 
[17,  (15)1, 
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= basic  transmission  loss  calculated  reference 
level  [17,  (17)], 


‘T 

Y(q) 


» a parameter  from  IF-73  [17,  (182)], 

= a long-term  variability  parameter  which  is 
calculated  using  (31)  and  (32)  with  c values 
from  tables  10  and  11, 

YeiCq)  = fgj,  Y(q).  (56) 

(O  if  lobing  option  [17,  sec.  3.1.1]  is 
,jused  and  the  aircraft  is  within  10  lobes! 
jof  its  radio  horizon 

r (L.. 


'YI 


(57) 


v^Hf  ■ ■ f’^br  ■ otherwise 


( 0 i f Ay  j ;l0  • 

*<^10  if  Ayj  ^10 

'Ayj  otherwise 


(58) 


Y (q<0.1)  = 


lesser  of 


lesser  of 


or  I for  lobing 


Yei(q) 

or 

^br'*'\  ■ ^Hf  ■ '^Y^j 


otherwise  . (59) 


Where  Cy  is  6,  5.3^  and  5 dB  for  q values  of  0.001,  0.001,  and 
0.01  respectively, 


' hf*  10  \ ^ 10 


YeCq*0*l)  “ (0.1)  otherwise 


and 


Yg(q  > 0.1)  = (60) 

These  equations  replace  similar  equations  in  IF-73  [17,  sec.  A. 5] 

- 7.  TUANvSITION  DISTANCB 

The  transition  distance  d^  is  used  in  blending  attenuation, 
valid  within  line-of-s.ight , with  the  radio  horizon  value.  It 
is  the  largest  distance  in  the  line-of-sight  region  at  which, 
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diffraction  effects  associated  with  terrain  are  considered  negli- 
gible. Values  estimated  for  d^  in  IF- 73  [17,  (140)]  have  been 
found  to  be  too  small  when  low  antennas  are  used  for  both  antennas. 
To  correct  this  difficulty,  d^  estimates  in  IF- 77  are  made  using 


dj^j^  when 

^x/6  ^h/6  "‘^Ll  h 

d,  otherwise 


(61) 


where  d^^^  is  the  horizon  distance  for  the  lower  terminal  (sec.  9.2), 


‘>/6 


is  the  distance  at  which  the  path  length  difference,  Ar 


[17,  (56)],  is  equal  to  X/6  (A  is  wave  length),  and  d^  is  the 
dp  of  IF-73  [17,  (140)].  The  distance  is  the  largest  distanc< 

at  which  a free- space  value  is  obtained  in  a two  ray  model  of 
reflection  from  a smooth  earth  with  a reflection  coefficient 

of  -1.  _ _ 

8.  FREE  SPACE  LOSS 

The  ray  length,  r,  term  of  the  free-space  loss  portion  of 
IF-73  [17,  (15)]  when  computed  via  the  IF-73  formulation 
[17,  r^  from  (54)  for  line-of-sight  or  path  distance  d for  trans- 
horizon paths]  can  give  loss  values  that  are  much  too  low  when 
a very  high  (satellite)  antenna  is  involved.  To  extend  IF-73 

to  such  cases  in  IF- 77,  a new  formualtion  for  r was  developed. 
This  formulation  may  be  summarized  as  follows: 

a^  = actual  earth  radius  (6370  km  = 3440  n mi) , 

d = great-circle  path  distance, 

^L1  distances, 

h. , 2 * horizon  elevations  (above  msl)  from  (72)  and 
* IF-73  (see  eqn.  396  of  App.  A). 

h^  2 “ antenna  elevations  (above  msl) , 

r^ » length  of  direct  ray  in  IF-73  [17,  (54)], 

r^H  = (h2-hj) ^+4 (h^+a^) (h2+ap) [sin(0 . 5d/a^) ] ^,  (62) 
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where  r^j^  is  the  within-the-horizon  ray  length  between  antennas 


above  an  air-less  earth  (i.e.,  bending  neglected)* 


'■U,2*  01.2-^11. (''•5di,i.2/^<,)]2(63) 


where  r^^j^  2 antenna  to  horizon  ray  length  for  an  airless 


earth  (no  ray  bending) , 


D. 


(64) 


where  D is  the  distance  between  horizons, 
s 


(65) 


where  rg^^  is  the  total  ray  length  for  a beyond-the-horizon  path, 


and 


r “ greater  of 


^o  ^WH 


d or  r 


BH, 


for  within-the-hor.izon 
or 

for  beyond-the-horizon 


(66) 


Equations  (62)  and  (63)  are  simply  the  application  of  the  half- 
angle law  of  cosines  formulation  where  two  sides  (a^^  + elevation) 
and  an  included  angle  (great-circle  distance  /a^)  are  known. 

9.  AIRBORNE  FACILITY 

The  IF- 77  version  allows  the  facility  (or  lower)  antenna  to 
be  airborne;  i.e.,  IF-73  was  extended  to  cover  air/air  and  air/ 
satellite  cases.  This  extension  involves  the  more  extensive  use 
of  parameters  based  on  ray  tracing  in  parts  of  the  model  associ- 
ated with  the  facility  antenna.  For  the  most  part,  these  para- 
meters are  similar  to  those  used  in  IF-73  for  the  aircraft  an- 
tenna only. 


9.1  Smooth  Earth  Horizons 

Ray  tracing  is  now  used  to  determine  the  smooth  earth  hori- 
zon distances  associated  with  both  terminals,  d^^^  2>  that  are 
used  in  the  calculations  associated  with  long-term  power  fading 
[17,  sec.  A. 5].  These  distances  are  determined  by  ray  tracing 
from  the  earth’s  surface  to  the  respective  antenna  heights  where 
the  initial  take-off  angle  is  0°  and  the  surface  refractivity 
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CN^  = 329  N-units)  corresponds  to  a 9000  km  (4860  n mi)  effective 
earth  radius. 

The  IF- 77  version  uses  ray  tracing  to  determine  smooth  earth 
horizon  distances  associated  with  both  terminals,  2* 

used  in  the  calculation  of  effective  antenna  heights.  These  dis- 
tances are  determined  by  ray  tracing  from  the  reflecting  surface 
elevation  [17,  fig.  13]  of  the  earth’s  surface  to  the  respective 
antenna  heights.  The  initial  take-off  angle  used  is  0°  and  the 


surface  refractivity , N^, 


is  calcualted  from  the  N or  effective 

0 


earth  value  specified  for  the  path  [17,  (18),  (20)].  Values  for 
dbsi  2 used  to  determine  effective  antenna  heights,  h^^  2» 
and  effective  heights  above  reflecting  plane,  Hj^  2»  follows: 

= effective  earth  radius  [17,  20]. 

= adjusted  earth  radius  [17,  (44)], 

= actual  earth  radius  (6370  km  = 3440  n mi) , 


a 

a 


‘al,2 


eg 


‘fc 


= actual  antenna  elevations  above  reflecting 
surface  elevation 

= height  of  facility  counterpoise  above  ground 
at  the  facility  site, 

■ height  of  the  facility  antenna  above  its  counter- 
poise, 

®sl,2  "*Lsl,2^“  (67) 


el,  2 


lesser  of 


“al,2 

or 


jo. 

L(a 


5 


[sec(eg2^  2^"^^  otherwi 


1,2  i “-I 

herwise 


Ah 


el, 2 


" ^al,2  ' ^el,2  » 


Ah 


al,2 


H, 


Ah 


h^j-  Ah^^for  ground  reflection 
h£^  for  counterpoise  reflection 


(68) 

(69) 

(70) 

(71) 
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”2  = 


ha2'^^a2^°^  ground  reflection 

h ,-Ah  -h  „ for  counterpoise  reflection 
a2  a2  eg 


These  expressions  are  extensions  of  similar  ones  used  in  IF- 73 
[17;  (34),  (32),  (45).  (46).  (48),  (49)]. 


9.2  Facility  Horizon 

The  IF- 73  version  allows  the  facility  horizon  to  be  specified 
by  (a)  any  two  horizon  parameters  (elevation,  elevation  angle,  or 
distance),  (b)  estimated  with  any  one  horizon  parameter  and  the 
terrain  parameter.  Ah,  (c)  estimated  from  Ah  alone,  or  (d)  calcu- 
lated for  smooth  earth  conditions  [17,  fig.  14],  Some  of  this 
flexibility  must  be  sacrificed  when  the  facility  is  high  since  the 
accurate  specification  of  more  than  one  horizon  parameter  requires 
prior  knowledge  of  ray  tracing  results. 

The  IF-  ' version  was  constructed  to  retain  all  facility  ho- 
rizon specification  flexibility  for  low  facility  antennas  and  yet 
allow  ray  tracing  to  be  used  for  high  facility  antennas.  This 
method  may  be  summarized  res  follows: 

1)  Determine  horizon  parameters  as  they  were  determined 
in  IF-73  [17,  fig.  14],  but  consider  the  results  as 
initial  values  that  may  be  changed  if  the  facility 
antenna  is  too  high.  The  resulting  parameters  are 


i m 


®Iei  ~ initial  estimate  of  the  horizon  elevation  angle 
®el’ 

hjLi  = initial  estimate  of  the  facility  horizon 
elevation  h^^^,  and 

*^111  ■ initial  estimate  of  facility  horizon  distance 
^ir 

2)  Facility  antenna  height,  h^,  and  effective  antenna 
height,  hg^,  from  (68)  are  used  to  tes:t  the  initial 


horizon  parameters  and  the  initial  parameter  values 
are  replaced  by  ones  appropriate  for  a smooth  earth 
if  the  test  conditions  are  met;  i.e.,  smooth  earth 
values  are  used  if 


hei  > 3 and 


Qt^i^O  and  h,>h 
lel  1 iLi 


>T«i<0  and  h,<h  \ 
lel-  1 iLi 


This  step  is  not  used  if  smooth  earth  parameters  were 
selected  in  step  2.  If  Ah^^^  from  (69)  is  zero  or 
less  the  initial  horizon  parameter  values  from  step  1 
are  used,  otherwise  ray  tracing  is  used  to  determine 


values  for  0^^^  and  dj^j^; 


•'Ll  ° ’'iLl 


®IeI  ''•'eli'’ 


Otherwise  use  ray  tracing 


‘^ILl  ^^el-° 


[otherwise  use  ray  tracing) 


The  ray  tracing  referred  to  in  the  equations  above  is  started  at 
the  horizon  elevation,  with  a take-off  angle  of  -0j^  and 

continues  until  the  facility  antenna  height  h^^  is  reached.  Then 
the  great-circle  distance  traversed  by  the  ray  is  taken  as 
and  the  negative  of  the  ray  arrival  angle  is  taken  as  0^^^.  The 
take-off  angle  used  is  calculated  from 


■®L  " '^®Iel'^^LlI^^^' 


10.  ANTENNA  PATTERNS 


This  section  deals  with  the  use  of  vertical  plane  antenna 
patterns  in  IF- 77.  These  patterns  give  gain  relative  to  the 
main  beam  gain  in  the  vertical  plane  [17,  sec.  A. 4. '2;  21,  figs. 
45  and  43].  Azimuth  patterns  are  used  only  in  program  TWIRL 
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(table  1)  which  is  described  in  the  APPLICATIONS  GUIDE  [21,  sec. 
3],  but  these  patterns  are  considered  to  be  more  a part  of  that 
particular  program  than  part  of  the  propagation  model  and  are  not 
discussed  here. 

10.1  Aircraft  Antenna 

The  aircraft  (or  higher  terminal)  antenna  pattern  in  IF- 73 
was  taken  as  isotropic,  and  modifications  to  allow  for  aircraft 
antenna  pattern  effects  are  included  in  IF- 77.  This  extension 
involves  the  use  of  the  gain  factors  which  are  discussed  in 
section  3.4. 

Aircraft  antenna  pattern  options  currently  built  into  the 
IF- 77  include  an  isotropic  antenna  and  a JTAC  [21,  (11)]  direc- 
tional pattern  where  the  half-power  beamwidth  and  the  tilt  of 
the  antenna  is  an  input  in  degrees.  Program  modifications  can 
easily  be  made  to  accommodate  other  patterns  that  are  specified 
in  terms  of  gain  versus  elevation  angle.  Horizontal  (or  azimuth) 
patterns  for  the  aircraft  antenna  are  not  used  in  any  of  the 
programs . 

Antenna  pattern  data  as  used  in  IF- 77  is  normalized  to  the 
main  beam  gain.  The  extent  to  which  the  main  beam  antenna  gain 
exceeds  that  of  an  isotropic  antenna  is  considered  as  a separate 
item  for  the  receiving  antenna  and  is  included  in  the  specifica- 
tion of  Equivalent  Isotropically  Radiated  Power  for  the  transmit- 
ting antenna  (EIRP) ; i.e., 

EIRP  [dBW]  = PYj^{dBW]  + G^[dBi]  (77) 

where  P.j.j^  is  the  total  power  radiated  from  the  antenna  and  is 
the  main  beam  gain  of  the  transmitting  antenna  relative  to 
isotropic. 

10.2  Ray  Elevation  Angles 

Incorporation  of  an  aircraft  antenna  pattern  into  the  model 
via  antenna  gain  factors  (sec.  3.4)  requires  that  gain  values  be 
obtained  from  the  pattern  for  the  direct  and  reflected  rays  at 
appropriate  elevation  Angles.  In  terms  of  the  variables  that 
are  similar  to  those  of  IF-73  [17,  sec.  A. 4. 2],  direct  ray 
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elevation  angles  0^^^^  2 ground  reflected  rays  2 given 
by 

(78) 

(79) 

(80) 


®hl,2=:^“"®l,2  ^ ®hP’ 

®H1,2  " ®hl,2  ^ ®L1,2, 

®gl,2  "®L1,2‘'^’®1,2» 


where  2 a smooth  horizon  elevation  angle  adjustment  term, 

and  the  remaining  parameters  are  calculated  as  in  IF-73;  i.e., 
o [17,  (53)],  0^  2 (50)T,  and  <|»  [77,  fig.  19].  Values-  for 

®L1  2 obtained  from 


®L1,2*  (®LR1,2  ■ ®LE1,2) 


(81) 


where  O^Rl  2 elevation  angles  of  the  smooth  earth  horizon 

rays  as  determined  with  ray  tracing,  2 elevation 

angles  of  the  horizon  rays  as  determined  using  the  effective 

earth  model,  and  the  remaining  parameters  are  calculated  as  in 

IF-73;  i.e.,  a^  [17,  (44)],  a [17,  (19)],  and  a [17,  (20)].  Values 
a o 

for  Gj^gj  2 are  obtained  from 

+a))  (82) 


®LE1,2"'^°^  ^^/^^1,2 


where  a [17,  (20)]  and  ^ (17,  (47)  (48)  i are  the  effective  earth 

radius  and  antenna  heights  of  IF-73.  ’Ine  effect  of  0^^^  2 

force  0u,  o and  0„i  to  have  the  values  obtained  via  ray  tracing 
nx I L gl % L 

at  the  smooth  earth  radio  horizon,  and  prorate  values  obtained 
elsewhere.  The  prorating  factor  (a^-a^) /(a-a^)  used  here  is  the 
same  factor  used  to  adjust  Ah^  [17,  (46)]  in  IF-73. 

10.3  Tracking  Options 

Tracking  options  are  available  for  each  terminal.  When  a 
tracking  option  is  used  for  a terminal,  its  antenna's  main  beam 
is  always  pointed  at  the  antenna  of  the  other  terminal  or  at 
the  radio  horizon  when  the  path  is  a transhorizon  path.  This 
is  accomplished  by  setting  the  beam  tilt  of  the  antenna  that  is 
tracking  to  the  direct  ray  elevation  angle  where  this  angle 
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becomes  the  horizon  elevation  angle  for  transhorizon  paths.  For 
example,  when  the  tracking  option  is  used  for  both  antennas  the 
full  gain  of  both  antennas  will  be  included  in  the  calculation 
of  transmission  loss  for  free-space  conditions. 

10.4  TACAN  Vertical  Pattern 

The  TACAN  RTA-2  vertical  pattern  used  with  IF- 77  [21,  fig. 

45]  is  based  on  a statistical  analysis  of  International  Telephone 
and  Telegraph  (ITT)  production  test  data  for  23  antennas.  Each 
of  these  antennas  were  tested  at  3 or  more  frequencies  so  that 
a total  of  78  patterns  were  used.  Of  these,  37  were  low  band 
(below  1088  MHz)  and  41  were  high  band  (above  1087  MHz) . Gain 
values  at  3“  intervals  for  elevation  angles  from  -60®  through 
60®  were  used  to  obtain  a distribution  of  gain  at  each  elevation 
angle.  Gain  values  exceeded  at  each  angle  for  5,  50,  and  95 
percent  of  the  data  are  shown  on  figure  9 along  with  the  standard 
deviation  of  the  gain  measurements.  Also  shown  in  figure  9 are 
gain  measurements  for  a single  antenna  in  the  above  60®  and 
below  -60®  range  that  were  obtained  from  a military  report  with 
limited  distribution,  and  the  piecewise  linear  approximation 
used  for  the  50  percent  in  IF-77. 

11.  SUMMARY 

This  report  covers  extensions  that  were  made  to  IF- 73  [17] 
in  the  process  of  developing  the  1977  capabilities  [21]  of  table 
1.  These  extensions  allow  the  programs  to  be  used  for  a wider 
variety  of  problems  such  as  those  involving  air/air  or  air/ 
satellite  propagation.  A brief  description  of  the  propagation 
model  provided  in  section  2 is  followed  by  detailed  discussions 
of  specific  model  extensions.  Minor  changes  to  and  errata  for 
IF- 73  are  provided  in  Appendix  A. 

The  1977  propagation  model  (IF-77)  has  been  incorporated  into 
computer  programs  that  are  useful  in  estimating  the  service  cover- 
age of  radio  systems  operating  in  the  frequency  band  from  0.1  to 
20  GHz.  They  may  be  used  to  obtain  a wide  variety  of  computer- 
generated microfilm  plots.  A plotting  capability  summary  is 


1;  -1^ 

I I-  - ■- 


TACAN  RTA-2 


provided  in  table  1 and  program  input  parameters  are  summarized 
in  tables  2 through  4.  These  tables  were  taken  from  an  APPLICA- 
TIONS GUIDE  [21]  for  the  programs  where  these  capabilities  and 
parameters  are  discussed  in  detail. 

Potential  users  should  1)  read  the  brief  description  of  the 
propagation  model  provided  in  section  2 to  see  if  the  model  is  ap- 
plicable to  his  problem,  2)  select  the  program(s)  whose  output(s) 
are  most  appropriate  from  the  information  given  in  table  1 [21, 
sec.  3],  3)  determine  values  for  the  input  parameters  given  in 
tables  2 through  4 [21,  sec.  4],  4)  request  a cost  estimate  for 
appropriate  computer  runs,,  and  5)  submit  the  formal  request  and/or 
purchase  order  that  may  be  required. 

Requests  to  the  FAA  should  be  addressed  to: 

Federal  Aviation  Administration 
Systems  Research  and  Development  Service 
Spectrum  Management  Staff,  ARD-60 
2100  Second  Street,  S.W. 

Washington,  D.  C.  20S91 

Attention:  Navigation  Specialist 

Telephone  contact  is  strongly  encouraged,  and  Mr.  Robert  Smith 
can  be  reached  at  426-3600  if  the  Federal  Telecommunications  Sys- 
tem (FTS)  is  used,  or  (202)  426-3600  if  commercial  telephone  is 
used. 

Other  requests  should  be  addressed  to: 

Department  of  Commerce 

Spectrum  Utilization  Division,  NTIA/ITS-1 
325  Broadway 
Boulder,  CO  80303 

Attention:  Mary  Ellen  Johnson 

Telephone  contact  is  strongly  encouraged,  and  Mrs.  Johnson  can  be 
reached  at  323-3587  if  FTS  is  used  or  (303)  499-1000  x 3587  if 
commercial  telephone  is  used.  If  extension  3587  can’t  be  reached, 
try  extension  4162,  which  is  the  Spectrum  Utilization  Division 
Office. 


APPENDIX  A.  CHANGES  FOR  FAA- RD- 73- 103 

Computer  Programs  for  Air/Ground  Propagation 
and  Interference  Analysis  0.1  to  20  GHz 


G.  D.  Gierhart  and  M.  E.  Johnson 


September  1973* 

All  changes  (errata  and  minor  modifications)  recommended  for  the 
above  report  by  the  authors  as  of  May  1978  are  listed  below. 

These  changes  ^ not  include  modifications  associated  with  the 
1977  extensions  that  are  discussed  in  the  text  of  the  present  re- 
port, but  do  include  some  minor  modifications  that  were  made  to 
accommodate  the  extensions.  Readers  finding  additional  errata  are 
urged  to  contact  an  author  at  the  U.  S.  Department  of  Commerce; 
Spectrum  Utilization  Division,  NTIA/ITS-1,  Boulder..  Colorado  80303. 
Mrs.  Johnson  can  be  reached  via  commercial  telephone  at  (303)  499- 
1000  X 3587  or  on  the  Federal  Telecommunications  System  (FTS)  at 
323-3587.  The  changes  are: 


Pa^e 


Location 


Line  2 


Line  1 


Changes 


End  of  first  . . . Change  the  (fig.  3)  to  (fig.  5). 
paragraph 


Change  the  Ns  to  . 

Change  "•..5o(q)  available  for  a fra..,- 

A 

tion  of  the  time  > q ...•’  to  "...S  (q) 
exceeded  for  a fraction  of  the  time 


*This  DOT  report  is  now  available  from  the  National  Technical 
Information  Service,  Operations  Division,  Springfield,  VA  22151. 
Order  using  accession  number  AD  770  335. 
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1 1% 


Pa'ge 


Location 


Chances 


Change  to 


©3  = 0.5(aVfD^''^/a  rad 
©4  = 1.5(a2/£)l/3/a  rad 


C90a) 


C90b) 


Fig.  20  caption  . 


(108) 

(118) 


(121) 


Change  the  last  h^^j^  to 

Change  the  dj^^  to 

Change  the  sign  of  h^^j^  2 minus 


“eLl,2 

Change  the  2.583  sin  (6^)  to 
5 . 1658  sin  (0.5  e^) . 


(126) Change  the  right  hand  side  to 


^ [(oeS^) V^L/^^^Ll^^KLZ^  . 


2 Sin 


(128) Change  the  right  hand  side  to 


Line  following 
(128)  . . . . 


(133)  . . . . 


-a  tan(0g)  +^(a  tan  0g)  + ?'^(l^Li‘^s2^ 


Insert  ’'hg2  from  (130)"  between 
and  "and." 

Change  the  2.583  sin  (0,)  to 
5.1658  sin  (0.5  6^).  ® 


(134)  Change  the  20  to  -20. 

(135)  Insert  "1"  between  "' 


and  "+". 


Second  line 

after  (135)  . 


(163) 


(182) 


Change  "...  e path..."  to  "...K 
path. . .". 

Change  the  right,  hand  side  to 
0.5696  hjj  { 1 + n exp  [-3.8  (0.1 

Change  the  right  hand  side  to 

L^ (0 . 5) - [L^£- Vg (0 . 5 , d^) - 20  log (R^g+Rtc^  ^ 
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iMl 

76 


86 


87 

104 


108 

111 


Location 
Line  2 
(184)  . . 


Line  5 


(194) 


(196) 


First  line 

after  (196) 


4th  line  after 
statement  22 


3rd  line  after 
statement  25  . 


Statement  8 . 
Line  5 . . . 


• • • * 


Changes 

Change  the  (IS)  to  (16). 

Change  the  L,  (0.5)  to  L,  +A 

br  y * 

Change  from  (66)"  to 

'’OhSin(,;,)/A  = 6-.^ 

Change  the  right  hand  side  to 

0.01  + 946  6^  if  6 < 0.00325 
6*15  6 if  0.00325  ^6<  0.0739 

0.45  +V^.000843-(6-0.10262  if  • . 
0.0739  <6<  0.1237 

0.601-1.06  6 if  0.1237  < 6 < 0.3 

0.01  + 0.875  exp  (-3.886)  otherwise 


Change  the  r^2  ^ r 2 r 2^  2 

S Q 


g 


■D 


Insert  ",  is  as  defined  for  (81)." 
between  "(40)"  and  "and  d". 


From  bottom 

lines  3 § 4 . 


"fOl'no'nf  (2)  to: 

'■rtl  I?”®*®’'?  specified"  and  ■ 

aie’speSfUd"!"®^'  elevation 


PMAY®®  for  PMIN, 

PMAX,  and  YC  to  "dB-W/sq  m".  . 


Change  the  U50  to  1)50. 


be  "HPFI"  and 

ISC"  should  be  "KE". 


Change  the  12  to  212, 


Change  "Read  8...IA"  to  "Read  8...IA, 
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^gaig^jatsaaHtaaa^^^ 


a??^)C5^5V.,  i 


m 


Page 


Location 


Changes 


Prom  page  bottom  Change  KE's  in  the  two  statements  pre- 
lines 7 § 8 ...  ceding  statement  73  (near  page  bottom) 

to  JD's. 


Immediately  after  Replace  the  whole  line  with  LE==11. 
first  comment 
statement 


After  statement 
136 


Insert  '’IFCSI.LE.SILIM.AND.ILB.IE.O.) 
GO  TO  137".  Attach  statement  numbers 
"138"  and  "139"  to  the  next  two  lines 
then  skip  a line  and  replace  the  next 
line  with  "IF(DZR.LT.O.)  GO  TO  145". 


Statement  46  . 


Replace  with  "46  RST= ( (RGP*RSP+RDG*RDG) / 
(GOD*GOD))+WA". 


Statement  after  Replace  with  "TLIM=+20 . *ALOG10 (GOD+REG+ 
statement  22  . REC) -GPD+(ADC18)*FTH)". 


After  statement  "137  WRL=CABS(G0D+AT2) 

135  insert  ....  IF(WRL.LE.GOD)  GO  TO  138 

WRL=CABS(G0D+AT1) 
WR*WRL*WRL+C.0901*GOD)$GO  TO  139. 


Right  before 
statement  148 


. Insert  145  DZR  = 0 . 


Line  7 Replace  ".195573364"  with  ".09679". 


Lin-;  before 

statement  30  . . . Delete 


Two  lines  after 

statement  30  . . . Replace  with  "TX= CCA*A/F)**THIRD)/A 

T3=0.5*TX  $ 1.5*TX". 


3rd  line  past 

statement  45  . . . Replace  with  "TH*2.*ASINF(CU*SQRTF 

CD/F*DL1*DL2)))". 


9th  line  after 

statement  45  . . Replace  "V5..."  with 

"V5=5.1658*SINFC0.5*TH5)*TMS". 


Line  6 Change  the  DLK4  to  D4. 


P^^Se  Location 


Changes 


Replace  *'V4=..."  with 

V4*S.1657(SINF(0.5*TH)*TM2)  .•' 

Statement  24  . . The  symbol  after  CALB  should  be  a "+»'. 

Replace  ’»PSWRB"  in  title  and  fi ret  n«o 
with  "PWSRB". 

Statement  24  . . insert  a "+*•  just  after  THET. 


APPENDIX  B.  LIST  OF  SYMBOLS 

This  list  includes  most  of  the  abbreviations,  acronyms,  and 
symbols  used  in  this  report  except  for  those  used  only  in  Appen- 
dix A.  Many  are  similar  to  those  previously  used  in  other  re- 
ports [17,  18,  21,  26,  35,  44].  The  units  given  for  symbols  in 
this  list  are  those  required  by  or  resulting  from  equations  as 
given  in  this  report.  Except  where  otherwise  indicated,  equations 
are  dimensionally  consistent  so  that  appropriate  units  can  be  se- 
lected by  the  user. 

In  the  following  list,  the  English  alphabet  precedes  the  Greek 
alphabet,  letters  precede  numbers,  and  lower-case  letters  precede 
upper-case  letters.  Miscellaneous  symbols  and  notations  are  gi- 
ven after  the  alphabetical  items. 

a Effective  earth  radius  as  calculated  in  IF- 73 

[17,  (20)]. 

app.  Appendix. 

An  adjusted  effective  earth  radius  shown  in 

figure  2 [17,  (44)]. 

p Principal  radii  of  curvature  of  reflecting  sur- 

face  of  the  reflecting  point  and  within,  a , or 
normal,  a , to  the  plane  of  incidence.  Usid 
in  (3).  " 

Sp  Actual  earth  radius  (6370  km = 3440  n mi). 

A A parameter  used  in  tropospheric  scatter  calcu- 

lations, from  (44). 

APODS  A program  name  (table  1) . 

ARD  Aviation  Research  and  Development. 

ATADU  A program  name  (table  1) . 

ATLAS  A program  name  (table  1). 

ATOA  A program  name  (table  1) . 

^pe  reference  plane  ray  bundle  area  (figure  1) 

associated  with  plane  earth  reflection  in  (2). 


C -'■‘S' 


A^Cq) 


A„(q) 


1,2,3 


'c,h,v 


Attenuation  [dB]  due  to  rain  calculated  via  (34) 
for  a fraction  of  time  q. 

Attenuation  [dB]  associated  with  rain  rate  and 
a fraction  of  time  q (sec.  4.4,  step  3). 

Terrain  attenuation  [dB]  from  (55)  that  is  as- 
sociated with  forward  scatter. 

The  reference  plane  ray  bundle  area  (figure  1) 
associated  with  plane  earth  reflection  in  (2). 

A conditional  adjustment  factor  [dB]  used  to 
prevent  available  signal  powers  from  exceeding 
levels  expected  for  free-space  propagation  by 
unrealistic  amounts,  from  (58). 

An  initial  of  value  of  A^  dB,  from  (57), 

Parameters  with  values  from  table  19  that  are 
used  in  (26). 

A parameter  used  in  tropospheric  scatter  cal- 
culation, from  (52) . i 


The  c-factor  from  tables  10  and  11  that  is  used 
in  (31)  and  (32). 


Centimeters  (10’ ^m). 

Phase  (rad)  of  plane  earth  reflection  coeffi- 
cient relative  to  ir  for  circular  (18) , horizon 
tal  (15),  and  vertical  (14)  polarization.  The 
total  phase  lag  associated  with  the  reflection 
coefficient  is  (tt-c  . ) rad. 

^ 9 A 9 V 

Parameters  with  values  from  table  19  that  are 
used  in  (26) . 

A parameter  used  in  tropospheric  scatter  cal- 
culations, from  (53). 


CCIR 


International  Radio  Consultative  Committee. 


CRPL 


Central  Radio  Propagation  Laboratory. 

Great  circle  distance  between  facility  and  air- 
craft. For  line-of-sight  paths,  it  is  calcula- 
ted as  indicated  in  figure  2. 

67 


Vi 


dB-W/sq  m 


Tsl,2 


'Ll,  2 


'Lo1,2 


DOC-BL 


DUDD 


Decibels,  10  log  (dimensionless  ratio  of 
powers . 

Antenna  gain  in  decibels  greater  than  isotropic. 

Power  in  decibels  greater  than  1 watt. 

Power  density  in  decibels  greater  than  1 watt 
per  square  meter. 

Degrees . 

The  dg  of  IF-73  [17,  (140)]  that  is  used  in  (61) 

to  calculate  d^  for  IF- 77. 
o 

Effective  distance  [17,  177]  that  is  used  in 
(26). 

An  initial  estimate  of  facility  horizon  dis- 
tance, made  via  IF-73  [17,  fig.  14]. 

Smooth  earth  horizon  distances  determined  via 
ray  tracing  (sec.  9.1). 

Horizon  distances  for  facility  and  aircraft 
respectively.  Values  for  d, , are  determined 
as  in  IF-73  [17,  (38)]. 

Smooth  earth  horizon  distances  determined  via 
ray  tracing  (sec.  9.1)  over  a 9000  km  (4860 
n mi)  earth. 

The  largest  distance  in  the  line-of-sight  region 
at  which  diffraction  effects  associated  with 
terrain  are  considered  negligible,  from  (61). 

The  largest  distance  at  which  a free-space 
value  of  basic  transmission  loss  is  obtained  in 
a two  ray  model  of  reflection  from  a smooth 
earth  with  an  effective  reflection  coefficient 
of  -1.  This  occurs  when  the  path  length  differ- 
ence, LX  [17,  (56)]  is  equal  to  X/6. 

Divergence  factor,  from  (4) . 

United  States  Department  of  Commerce,  Boulder 
Laboratories  ^ 

United  States  D^bj/artment  of  Transportation. 

A program  name  (table  1) . 


DURATA  A program  name  (table  1) . 

D/U  Desired- to-undesired  signal  ratio  [dB]  avail- 
able at  the  terminals  of  an  ideal  (loss  less) 

receiving  antenna. 

Distance  between  radio  horizons,  from  (64). 

D-  , Distance  shown  in  figure  2 [17,  (51)]. 

eqn.  Equation. 

exp(...)  Exponential;  e.g.,  exp(2)  « e^ 

EIRP  Equivalent  Isotropically  Radiated  Power 

TdBW]  " ~ " 

ESSA  Environmental  Science  Services  Administration, 

f Frequency. 

fss  Facility  £ite  £urface  (table  2) , 

ft  Feet. 

ft- fss  Feet  above  facility  site  surface, 

ft-msl  Feet  above  mean  sea  level. 


FAA 

FTS 

^dcfh 


g 

g(g,f) 

®D,R 


Elevation  angle  correction  factor  [17,  (179)]. 

Federal  Aviation  Administration. 

Federal  T^elecommunications  System 

Reflection  reduction  factor  associated  with 
diffuse  reflection  and  surface  roughness 
[17,  (194)]. 

Reflection  reduction  factor  associated  with 
ray  lengths,  from  (7). 

Specular  reflection  reduction  factor  associated 
with  surface  roughness  [17,  (66)]. 

Normalized  voltage  antenna  gain  used  for  the 
facility  antenna  in  IF- 75  [17,  (67)]. 

Frequency  gain  factor  from  (29)  or  (30) . 

Voltage  gain  [V/V]  factors  associated  with 
direct  and  reflected  rays,  from  (10)  and  (11). 
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Sd1,2 


®hDl,2 

®hRl,2 


'Rl,2 


®vDl,2 
®vRl , 2 


GOES 

Gn-i  / 


"Rl,2 


^al,2 

If. 

^cg 

4/^  r ■ 

^el,2 

' \ 

>>.0 

f:fS' 

^ILl 

Voltage  gain  [V/V]  of  terminal  antennas  in  the 
direction  of  the  direct  ray  (figure  3)  relative 
to  main  beam  gain. 


Voltage  gain  [V/V]  similar  to  gp,  2» 
specifically  for  horizontal  polarization. 


Voltage  gain  [V/V]  similar  to  gj,,  ,» 
specifically  for  horizontal  * 


polarization. 


Gain  factor  for  the  reflected  ray  and  vertical 
polarization,  from  (12). 


Gain  factor  for  the  reflected  ray  and  horizon-' 
tal  polarization,  from  (13). 


Voltage  gain  [V/V]  of  terminal  antennas  in-  the 
direction  of  the  reflected  ray  (figure  3)  rela- 
tive to  main  beam  gain. 


Voltage  gain  [V/V]  similar  to  "g^,  speci- 

fically for  vertical  polarization: 

Voltage  gain  [V/V]  similar  to  gj.,  - speci- 
fically for  vertical  polarization: 


Gigahertz  (10^  Hz) 

Geostationary  Operational  Environmental  Satellite. 


Gains  [dB]  g„,  ..  expressed  in  decibels,  from 
(9). 


Main  bieam  gain  [dBi]  of  transmitting  antenna. 
Hours . 


Actual  antenna  elevations  above  reflecting  sur- 
face elevation. 


Height  of  facility  counterpoise  above  ground 
at  the  facility  site. 


Effective  antenna  height  above  h^^,  from  (68). 


Height  of  the  facility  antenna  above  its 
counterpoise. 


Initial  value  of  h. 


f-  -^< 


'Ll,  2 


“1,2 

HDBK 


HIPOD 


IEEE 


IF-73 


IF-77 


JTAC 


1,2 

LOBING 


Terminal  horizon  elevations. 


Elevation  o£  reflecting  surface  above  msl. 


Common  volume  elevation  above  average  terrain. 


Antenna  elevations  above  msl, 


Handbook. 


A program  name  (table  1) . 


Hertz. 


Significant  wave  height  (table  5) 


Antenna  elevations  shown  in  figure  2 [17,  (52)]. 


An  index  for  specific  transmission  loss  levels 
used  in  the  distribution  mixing  process  (sec. 
4.1).  It  has  values  from  1 to  M. 


Inches . 


institute  of  Electrical  and  Electronic  Engineers. 


Institute  for  Telecommunication  Sciences. 


institute  for  .Radio  Engineers, 
international  Telephone  and  Telegraph. 

ITS-FAA-1973  propagation  model. 
ITS-FAA-1977  propagation  model. 


/-I  or  an  index  (1  to  N)  for  specific  trans- 
mission loss  distributions  used  in  the  dis- 
tribution mixing  process  (sec.  4.1). 


Joint  Technical  Advisory  Committee. 
Kilometer  (10 ^m). 


Total  ray  length,  from  (45). 


Pommon  (base  10)  logarithm. 


Terminal  to  common  volume  ray  lengths. 


A program  name  (table  1) , 


} 


mhos 


n mi 


nsec 


NOAA 


NTIA 


NTIS 


N-unit 


Basic  transmission  loss  [dB]  for  free  space 
[17,  (15)]. 

Basic  transmission  loss  [dB]  calculated  refer- 
ence level  [17,  (17)]. 

Meters . 

Unit  of  conductance  or  siemens. 

Minutes . 

Millimeters  (10‘^m). 

Mean  sea  level. 

Number  of  transmission  loss  levels  used  in 
mixing  distributions  which  is  also  the  final 
value  for  the  index  i (sec.  4.1). 

Military . 

Megahertz  (10*^  Hz), 

A power  used  in  the  frequency  scaling  via  (37) . 
Nautical  miles. 

9 

Nanoseconds  (10"  sec) . 

Number  of  distributions  to  be  mixed  which  is 
also  the  final  value  for  the  index  j (sec.  4.1), 
or  North  latitude. 

National  Bureau  of  Standards. 

National  Oceanic  and  Atmospheric  Administration. 

National  Telecommunications  and  Information 
Administration. 

National  Technical  information  S^ervice. 

Minimum  monthly  mean  surface  refractivity 
(N- units)  referred  to  mean  sea  level  [17, 
figure  3] . 

Minimum  monthly  surface  refractivity  in  N-units 
[17,  (18)]. 

Units  of  refractivity  [4,  sec.  1.3]  correspond- 
ing to  10®  (ref  active  index  -1). 


I 
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TR 


‘^1.2 


rad 

rms 

^BH 


eo,s,w 


LI, 2 


"WH 


1,2 


12 


Total  power  [dBW]  radiated,  used  in  (77). 

Dimensionless  fraction  of  time  used  in  time 
availability  specification;  e.g.  Yq  (0.1) 
where  q = 0.1  implies  a time  availability  of 
10  percent. 

Parameters  used  in  tropospheric  scatter  calcu- 
lations from  (51) . 

Time  availabilities  for  mixed  distributions 
that  correspond  to  specific  transmission  loss 
levels,  from  (22). 

Time  availabilities  for  each  transmission  loss 
level  (index  i)  of  each  transmission  loss  dis- 
tribution (index  j)  involved  in  the  distribu- 
tion mixing  process  (sec.  4,1). 

Ray  length  used  in  the  calculation  of  free 
space  loss,  from  (66). 

Radians . 

Root  mean  square. 

Ray  length  for  beyond -the -horizon  paths,  from 
(65). 

Effective  ray  lengths  (sec.  4.4)  for  attenua- 
tion associated  with  oxygen  absorption,  (r^o) » 
rain  storm  attenuation  (r^s) , and  water  vapor 
cbsorption  (r^^) . 

Direct  Ray  length  shown  in  figure  2 [17,  (54)]. 

Antenna  to  horizon  ray  lengths  for  airless 
earth,  from  (63). 

In- storm  ray  length  used  in  rain  attenuation 
calculation,  from  (33) , 

Within- the-horizon  ray  length  for  airless  earth, 
from  (62). 

Segments  of  reflected  ray  path  shown  in  figure 
2 and  components  of 

Reflected  ray  path  length  as  shown  in  figure  2 
[17,  (55)]. 


73 


RTA-2 


c,h,v 


Magnitude  of  complex  plane  earth  reflection 
coefficient. 

A TACAN  antenna  type. 

Magnitudes  of  complex  plane  earth  reflection 
coefficients  for  circular,  horizontal,  and 
vertical  polarization. 

A parameter  used  in  the  calculation  of  the  di- 
vergence factor,  from  (5). 

Modules  of  asymmetry  used  in  tropospheric  scat- 
ter calculations,  from  (43). 

Seconds . 


sq  m 


S Jal 


Square  meters. 
Statute  miles. 


SRVLUM 


Super  High  Frequency  (3  to  30  GHz) . 
Inverse  sine  with  principal  value, 

A program  name  (table  1) . 


TACAN 


TWIRL 


eo,s,w 


Scattering  efficiency  term  [dB]  used  in  tropo- 
spheric scatter  calculations,  from  (42). 

Scattering  volume  term  [dB]  of  tropospheric 
scatter  calculations,  from  (54). 

Relaxation  time  [ys]  used  in  the  calculation 
of  surface  constants  for  water  (table  6). 

TACtical  Air  Navigatio...  an  air  navigation  aid 
used  to  provide  aircraft  with  distance  and 
bearing  information. 

A program  name  (table  1) , 

Height  or  layer  thickness  (sec.  4.4)  used  in 
attenuation  calculations  for  oxygen  absorption 
(T  ),  rain  storm  attenuation  (T  ) , or  water 
vapor  absorption  (T  ) . 

Ultra-High  Frequency  (300  to  3000  MHz) . 

Volts. 


I Pc. 

* ' //*  /• 


VJq) 

V 

l.i.M 

W 

Wi  s M 

YCq) 


i 

! 


i 


if,(q) 

’fel<‘>) 

VjCq) 

Yic(q) 

Vj(q) 

^T 

Yo(O.l) 

Y„C0.9) 

’^136^''' 

Y,(q) 

Y,(q) 


Variability  for  specific  climate  or  time  block, 
from  (23). 


Variability  levels  (Vj , . . . Vi , . . . V^)  used  in 
mixing  process  (sec.  4.1). 

Watts . 


Weighting  factors  (Wj , . . .Wj , . . . W^)  used  in 
mixing  process  (sec.  4.1). 

Variability  (dB  greater  than  median)  of  hourly 
median  received  power  about  its  median,  from 
(27),  (28),  (31)  and  (32). 

A complex  parameter  used  in  the  calculation  of 
the  plane  earth  reflection  coefficient,  from 
(16). 

Effective  variability  (dB  greater  than  median) 
of  hourly  median  received  power  about  its 
median,  from  (59)  and  (60) . 

Initial  value  Y (q)  from  (56). 

0 

Variability  [dB]  associated  with  ionospheric 
scintillation  (figure  7) . 

y,(q)  for  a particular  distribution  to  be  used 
in  the  mixing  process  to  obtain  resultant  Yj(q), 
from  (25). 

Variability  (dB  greater  than  median)  associated 
with  rain  attenuation,  from  (35). 

A parameter  from  IF-73  [17,  (182)]. 

A reference  variability  level  used  to  calculate 
Y(O.l),  from  (26). 

A reference  variab'lity  level  used  to  calculate 
Y(0.9),  from  (26). 

Variability  associated  with  ionospheric  scin- 
tillation at  136  MHz. 

Variability  (dB  greater  than  median)  associated 
with  multipath  [17,  p.  38], 

Total  variability  (dB  greater  than  median),  from 
(36). 


An  angle  shown  in  figure  2 [17,  (53)]. 

A parameter  in  tropospheric  scatter  calculations, 
from  (40) . 

Terrain  parameter  used  to  characterize  terrain 
[17,  sec.  A. 4.1;  26,  sec.  2.2]. 

Adjusted  effective  altitude  correction  factors, 
from  (70). 

Effective  altitude  correction  factors,  from 
(69). 

Path  length  difference  for  rays  shown  in  figure 
2 [17,  (56)]. 

Dielectric  constant  from  table  6 or  calculated 
for  water  using  (20). 

Complex  dielectric  constant,  from  (17), 

Dielectric  constant  representing  the  sum  of 
electronic  and  atomic  polarizations.  For 
water,  = 4.9, 

Static  dielectric  constant  (table  6). 

Parameters  used  in  tropospheric  scatter  calcu- 
lations, from  (39)  and  (41). 

A parameter  used  in  tropospheric  scatter  calcu- 
lations, from  (46). 

Scattering  angle  used  in  tropospheric  scatter 
calculations.  It  is  the  angle  between  trans- 
mitter horizon  to  common  volume  ray  and  the 
common  volume  to  receiver  horizon  ray  as  both 
leave  their  crossover  point. 

Elevation  angle  of  horizon  from  the  facility, 
from  (74). 

Elevation  angle  of  horizon  from  the  aircraft 
[17,  (39)]. 

Latitude  of  earth  facility  for  (38) . 

Elevation  angles  of  the  ground  reflected  rays 
at  the  terminal  anteunrs,  from  (80). 

Parameter  used  to  calculate  9jjj^  2 f^'om  (78). 

76 


0 


lel 


nirnct  ray  elevation  angles  at  the  terminal 
antennas,  from  (79). 

Initial  estimate  of  0 , ; i.e.,  0^,  as  calcula- 
ted in  IF-73  [17,  figfife  14]. 


0 


L 


A ray  tracing  take-off  angle  at  the  facility 
horizon,  from  (76). 


0 


LI, 2 


Horizon  elevation  angle  adjustment  terms, 
from  (81) . 


®LE1,2 


Horizon  elevation  angles  as  determined  with  ? 

effective  earth  radius  model,  for  (81).  I 


®LR1,2 

®sl,2 


Horizon  elevation  angles  as  determined  with 
ray  tracing,  for  (81) . 

Central  angles  below  the  smooth  earth  terminal  to 
horizon  distances  for  the  effective  earth 
model,  from  (67). 

Angle  defined  in  figure  2. 


®1,2 

K 

A 

liS 

^1.2 


^1,2 


Angles  shown  in  figure  2.  I 

Wave  number,  from  (49).  I 

I 

Wavelength.  ^ 

Microseconds  (10'®  sec).  | 

Parameters  used  in  tropospheric  scatter  ca'lcu-  | 

lations,  from  (50).  I 

Surface  conductivity  [mho/m]  from  (21)  or  ? 

table  6.  ! ] 


Root-mean- square  deviation  of  surface  excur-  | 
sions  within  the  limits  of  the  first  Fresnel  | 
zone  in  the  dominant  reflecting  plane  [17,  | 
(65)],  from  table  5 or  (1).  I 

5 


Ionic  conductivity  [mho/m],  from  table  6. 


Parameters  used  in  tropospheric  scatter  calcu- 
lations, from  (47)  and  (48). 

Grazing  angle  shown  in  figures  2 and  3. 

Grazing  angle  associated  with  the  pseudo  Brew- 
st-  r angle,  from  (19). 


I 

I 


(...) 


Degrees  Celsius. 
Degrees  fahrenheit. 
Degrees;  e.g.  12®. 


Expression  evaluated  for  specific  conditions 
such  as  climate  or  time  block  in  (23) . 

Expression  evaluated  for  radio  horizon  conditions, 
Magnitude  of  expression;  e.g.,  ll"jl»/2~ 
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